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ABSTRACT 

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) com- 
prises an Hl-selected sample of 155 galaxies spanning all morphological types with 
distances less than 25 Mpc. We describe the scientific goals of the survey, the sample 
selection, and the observing strategy. We also present an atlas and analysis of the 
CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer 
Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions 
in the galaxies in this sample and explore the correlation of the far-infrared luminosity, 
which traces star formation, with the CO luminosity, which traces the molecular gas 
mass. By comparing the NGLS data with merging galaxies at low and high redshift 
which have also been observed in the CO J=3-2 line, we show that the correlation of 
far-infrared and CO luminosity shows a significant trend with luminosity. This trend 
is consistent with a molecular gas depletion time which is more than an order of mag- 
nitude faster in the merger galaxies than in nearby normal galaxies. We also find a 
strong correlation of the Lfir/ Lco(3-2) ratio with the atomic to molecular gas mass 
ratio. This correlation suggests that some of the far-infrared emission originates from 
dust associated with atomic gas and that its contribution is particularly important in 
galaxies where most of the gas is in the atomic phase. 

Key words: galaxies: ISM — galaxies: kinematics and dynamics — galaxies: spiral 
- ISM: molecules - stars: formation 



1 INTRODUCTION 

Star formation is one of the most important processes driv- 
ing the evolution of galaxies. The presence or absence of 
significant star formation is one of the key characteristics 
which distinguishes spiral and elliptical galaxies. The intense 
bursts of star formation triggered by galaxy interactions and 
mergers produ ce some of the most lumin ous galaxies in the 
local universe (|Sanders fc Mir abel 199f|). At high redshift, 
many galaxies are seen to be forming stars at rates which 
far exceed those of all but th e most extreme local merg- 
ers (|Tacconi et alj [20081 . 12009T ). Since stars form from gas, 
specifically from molecular gas, understanding the proper- 
ties of the interstellar medium (ISM) is critical to under- 
standing the rate and regulation of st ar formation in galaxies 
l|Lerov et al.ll2008l : IBigiel et alj|2008h . 

At the most basic level, the amount of gas in a galaxy is 
an important constr aint on the amount of star formation the 
galaxy can sustain dKennicut t 1989; Kennicutt et al. I l2007l : 



iLerov et al] 120081 : IBigiel et al 



20081 ) . Recent studies have 



shown that it is the amount of molecular gas that is most 
impo r tant, rather than the total gas cont ent (|Bigiel et al.l 
|200S| ; ILerov et al]|200Sl ; IBigiel et al]|201ll ). This picture is 
consistent with Galactic studies which show that stars form 
exclusively in molecular clouds, and most commonly in 
the densest regions of t hose clouds l|Lada et al.1 Il991al lbl: 
lArzoumanian et alj|201ll ). This suggests that the proper- 
ties of the molecular gas, in particular its average density 
and perhaps the fraction of gas above a critical density on 
sub-parsec scales, are likely to affect the resulting star for- 
mation. Star formation in turn can affect the properties of 
the dense gas, by increasing its temperatu re (| Wilson et al.l 
ll997l : lMeier et al.ll200ll ; iTosaki et al.ll2007l ) a nd perhaps by 
trigg ering a subsequent generation of stars (jZavagno et al.l 
|2010| ). Unusual environments are also likely to affect both 
the gas properties and the star formation process. High shear 
in galactic bars, harassment in a dense galaxy cluster, and 
galaxy mergers and interactions have the potential to either 



dampen or enhance the star formation process. The wide 
range of environmental processes at work, both on galac- 
tic and extragalactic scales, implies that large samples of 
galaxies are required to tease out the most important ef- 
fects, while high resolution is required to isolate individual 
star forming regions, separate arm from interarm regions, 
and resolve galactic bars. 

There have been a number of large surveys of the atomic 
gas content of nearby galaxies, of which some of the most 
recent include The Hi Nearby Galaxy Survey, THINGS 
l| Walter et al.ll2008l). the Arecibo L egacy Fast ALFA survey, 
ALFALFA dGiovanelli et al.ll2005t) . and the VLA Imaging 
of Vi rgo Spirals in Atomic Gas survey, VIVA l|Chung et al.l 
l2009t ). However, only the 34 galaxies in the THINGS and 
the 53 galaxies in the VIVA sample have sufficient spatial 
resolution to probe scales of one kiloparsec and below. Com- 
pared to the Hi 21 cm line, the CO lines used to trace 
molecular gas are relatively more difficult to observe due 
to their shorter wavelengths and the smaller field of view 
of millimeter-wave radio telescopes equipped with single 
pixel detectors. As a result, most CO extragalactic surveys 
have sa mpled a relatively small region at the c enter of each 
galaxy ilYoung et al.lll995l : iHelfer et all 120031 : iBraine et af 
I1993I ; iDumke et al.ll200ll ). Two recent surveys (TKuno et al 



2007); ILerov et al.f 2009 ) have used array receivers to observe 
large unbiased regions, although still in relatively small (18 
and 40) samples of galaxies. Finally, dust continuum obser- 
vations in the submillimeter to far-infrared present an alter- 
native method of tracing the molecular gas that requires a 
mostly independent set of physical parameters, such as dust 
emissivity, gas-to-dust mass r atio, and the atomic gas mass 
ilThronsonlll98Sl ; Ilsraellll997l : lEales et al]|2010l : ILerov et al.1 
120111 ). Two large surveys have been made at 850 /im, one of 
an Infrared Astron omical Satellite (IRAS) selected sample 
Dunne et allli oOO) and one of an optically selected sample 



Vlahakis et al.ll2005l ). The galaxies selected for these 



pies were relatively distant (> 25 Mpc) and thus primarily 
global measurements of the dust luminosity were obtained. 
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The Herschel Reference Survey (|Boselli et al.ll2O10l ) is ob- 
serving 323 galaxies with distances between 15 and 25 Mpc 
at 250, 350, and 500 /im and has significant overlap with 
the sample presented here. Other Herschel surveys which 
overlap with the NGLS sample include the Key Insights 
on Nearby Galaxies: A Far-Infrared Survey with Herschel 



jKennicutt et alj 



20 111) and the Herschel Virgo Cluster Sur- 



vey (|Davies et al 



20ld ). Using the dust continuum to mea- 



sure the star forming gas is a promising avenue to explore, 
especially since this method is one that can in principle be 
used for galaxies at higher redshifts. 

Taking advantage of new instrumentation for both spec- 
tral line and continuum data on the James Clerk Maxwell 
Telescope (JCMT), we are carrying out the Nearby Galax- 
ies Legacy Survey (NGLS)Q, a large survey of 155 nearby 
galaxies using the CO J=3-2 line and continuum observa- 
tions at 850 and 450 fim. The survey is designed to address 
four broad scientific goals. 

(i) Physical properties of dust in galaxies. The continuum 
data from this survey will probe a range of the dust spec- 
tral energy distribution (SED) that is critical to determin- 
ing the total mass of dust as well as the relative propor- 
tion and physical properties of the different dust components 
(polycyclic aromatic hydrocarbons, very small grains, large 
grains). Most importantly, these long wavelength data will 
trace any excess submillimetre emission that causes an up- 
turn in the dust SED. The submillimetre ex cess could orig- 
inate from : dust with <10 K temperatures (|Galliano et al.l 
l2003l , l2005l .l 201 ll ); very small grains with shallow dust emis- 
sivities that a re not prominent in the dust SED betw een 60 
and 500 /im (jLisenfeld et alJ 120021 : IZhu et all h00$ ): large 
dust grains with enhanced dust emissivi ties at >850 /im 
jBendo et alj|2006l: lO'Halloran et al ]|201Ch : or spinning dust 
grains (|Bot et al.ll2010l : lAde et al.ll201ll ). 

(ii) Molecular gas properties and the gas to dust ratio. 
Our high-resolution data will allow us to compare the radial 
profiles of the dust, Hi, and CO emission within our well- 
selected sample. The CO J=3-2 line will effectively trace 
the warmer, denser molec ular gas that is more directly in- 
volved in star formation l|lono et al.ll2009l ). In comparison, 
the CO J=l-0 line also traces more diffuse and low density 
gas (|Wilson fc Walkerlll994l : lltosolowskv et al.ll2007l ). With 
observations of all three components of the ISM (molecular 
gas, atomic gas, and dust), we will be able to determine 
accurate gas-to-dust mass ratios and provide constraints 
on the variat i on of the CO-to-H2 conver sion factor Xgo 
JStrong et all Il988h bv fittin g the data (jThronsonl 1 19881 : 



llsraellll997l : lBraine et atlll997f ). 



(iii) The effect of galaxy morphology. The larger-scale 
galaxy environment can play a significant role in the proper- 
ties and structure of the dense ISM. For example, an increase 
in the density in the ISM in the centers of spiral galaxies 
has b een attributed to an increased pressure (|Helfer fc Blitd 
1993). Elliptical galaxi es have relatively l o w column den- 
sities of gas and dust l|Knapp et al.l 1 19891 : iBregman et al.l 
1988) combined with an intense radiation field dominated 
by older, low-mass stars and, in some cases, substantial X- 
ray halos. Early-type galaxies also show more compact and 
symmetric distributions of dust at < 70 fim compared to 



http: / /www .jach.hawaii.edu/JCMT/surveys/ 



late- t ype galaxies (|Bendo et al.l 120071 : iMunoz-Mateos et al.l 
2009). In spiral galaxies, dust and gas p roperties may dif- 
fer between arm and inter-arm regions (jAlton et all I2OO2I : 
iFovle et al.ll201Cl ), while the role of spiral arms in the star 
formation process is the subject of considerable debate 
llElmegreen fc Elmegreenlll986l : IVogel. Kulkarni. fc Scovillel 
119881 ). With observations across the full range of galaxy mor- 
phologies and with sub-kiloparsec resolution, we will be able 
to study the effect of morphology on the molecular gas prop- 
erties. 

(iv) The impact of unusual environments. Metallicity has 
been shown to aff ect the structural properties of the ISM 
ijLerov et al.ll201ll ). Lower self-shielding is expected to pro- 
duce smaller regions of cold, dense gas that can be traced 
by CO emission wi th relatively larger and wa rmer photon- 
dominated regions ( Mad den et al 1 1 1991 1 20061 ). Galaxies re- 
siding in rich clusters ca n be affected by ram pressure 
stripping dLee et al l 120031 ) and gravitational harassment 
IIMoore et al.l 1 19981 ) which reduces their ISM content rela- 
tive to field galaxies. The halos of spiral galaxies can be 
populated with gas via superwinds or tidal interactions 
ijRoussel et al.ll20ld) or even by r elatively normal rates of 
star formation |Lee fc rrwinl fl997). With its galaxies span- 
ning the full range of environment, from isolated galaxies to 
small groups to the dense environment of the Virgo cluster, 
the NGLS will be able to quantify the effect of environment 
on the molecular ISM in galaxies. 

Most surveys of molecular gas in our own or other 
nearby galaxies have used the ground state CO J=l-0 line 
as a tracer. The choice of the CO J=3-2 line for the NGLS 
was driven by the available instrumentation at the JCMT. 
Compared to the CO J=l-0 line, (5.5 K above ground with 
a critical density of 1.1 x 10 3 cm" 3 ), the CO J=3-2 line (33 
K and 2.1 x 10 4 cm -3 ) traces relatively warmer and denser 
gas. (Since both lines are usually optically thick, the effective 
critical density is likely reduced by a factor of 10 or more.) 
There is growing evidence that the CO J=3-2 emission cor- 
relates more tightly with the star formation ra te or star for- 
mation efficiency than doe s the CO J=l-0 l ine (|Wilson et al.l 
l2009l : lMuraoka et al.ll2007l ). iKomugil l|2007f ) showed that the 
CO J=3-2 emission correlates linearly with star formation 
rate derived from extinction-corrected Ha emission and that 
the correlation was tight er with th e J=3- 2 line than with 
the J=l-0 line. Similarly, llono et all (|2009l ) showed that the 
CO J=3-2 emission correlates nearly linearly with the far- 
infrared luminosity for a sample of local luminous infrared 
galaxies and high redshift submillimeter galaxies. Thus, it is 
appears that the CO J=3-2 emission is preferentially tracing 
the molecular gas associated directly with star formation, 
such as high-density gas that is forming stars or warm gas 
heated by star formation, rather than the total molecular 
gas content of a galaxy. 

In this paper, we describe the NGLS sample selection, 
the CO J=3-2 observations and data reduction, and the 
planned observing strategy for continuum 450 and 850 fim 
observations (J2). In § [3] we present the CO J=3-2 inte- 
grated intensity images, as well as maps of the velocity field 
and velocity dispersion for those galaxies with sufficiently 
strong signal. In § [4] we examine the molecular and atomic 
gas masses for the Spitzer Infrared Near by Galaxies Survey 
(SINGS) sample jKennicutt et al.ll2003l ) and compare the 
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CO J=3-2 luminosity with the far-infrared luminosity for 
both the galaxies in the SINGS sample, local luminous and 
ultraluminous infrared galaxies, and high redshift quasars 
and submillimeter galaxies which have also been observed 
in the CO J=3-2 transition |lono et al.ll2009l ). We give our 
conclusions in Previous papers in this seri es have exam- 



ined individual or small samples of ga laxies (Wilson et al 



20091; IWarren et al.ll2010l; iBendo et al.|[2010a; 



Wilson et al 



201ll ; llrwin et al.ll201ll ; ISanchez-Gallego et al.ll201ll ) Future 
papers will exploit the spatially resolved nature of these data 
by examining the CO line ratios and excitation (Rosolowsky 
et al. in prep.) and the gas depletion times (Sinukoff et al., 
in prep.). 



2 OBSERVATIONS AND DATA PROCESSING 
2.1 Sample selection 

The galaxies in the NGLS are an Hi-flux selected sample. 
This selection method was chosen to avoid biasing the sur- 
vey towards galaxies with higher star formation rates, as 
might be the case for far-infrared or blue magnitude selec- 
tion criteria, while still targeting galaxies with a significant 
interstellar medium, as m ight not be the cas e for a purely 
mass selected sample, e.g. iBoselli et al] (|2010t ). Our sample 
should also not be unduly biased by dust content being ei- 
ther particularly hot or cold, and should be reasonably rep- 
resentative of the local submillimetre universe. H I emission 
typically extends to much a larger radius than CO or stellar 
emission, so that our sample may include galaxies that are 
Hi-rich but HVpoor. 

To obtain good spatial resolution while still being able 
to obtain sensitive maps in a reasonable period of time, 
the galaxies were selected to have distances between 2 and 
25 Mpc. This distance limit excludes galaxies in the Local 
Group but does include galaxies in the Virgo cluster, which 
allows us to test the effect of cluster environment on the 
dense ISM. We used the HyperLedc0 data base in February 
2005 to extract all non- Virgo galaxies (see below) with H I 
fluxes > 3.3 Jy km s _1 and Virgo-corrected galactocentric 
velocities v < 1875 km s . Out to these limits, the Hyper- 
Leda database is essentially complete. We removed galaxies 
with declinations below -25°; we further removed galaxies 
with Galactic latitudes between -25° and 25° to minimize 
the impact of Galactic cirrus in complementary data from 
Spitzer and Herschel. Virgo cluster member galaxies were 
selected to be galaxies within an 8° x 16° ellipse centered on 
M87 (RA=12.4h, DEC=12. 4°) and with v elocit ies between 
500 and 2500 km s" 1 , e.g. iDavies et all l|2004l) : all Virgo 
galaxies with an H I flux entry were included in this initial 
selection. These selection criteria gave us a sample of 1002 
field galaxies and 148 Virgo cluster galaxies. 

Since we estimated that at least 1500 hours of telescope 
time would be required to observe this complete sample of 
1150 galaxies, we needed to apply further s election criteria. 
We f irst selected all the SINGS galaxies (|Kennicutt et al] 
120031 ) that met our declination and Galactic latitude cri- 
teria; the wealth of complementary data avail able on these 
galaxies from the SINGS and other surveys l|Walter et al.l 
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2008) make them natural targets. Forty-seven SINGS galax- 
ies met our criteria (Tables [1] and [2j . Two SINGS galaxies, 
NGC 5194 (M51) and NGC 4789A (DDO 154), were ob- 
served in separate programs and were not re-observed as 
part of our survey; we present maps made from those data 
in this paper. 

We then estimated the minimum size of the field and 
Virgo cluster samples that would allow us to achieve the sci- 
entific goals of this study. For statistical studies of galaxy 
properties, we want to divide the galaxies into 4 morpho- 
logical bins (E/S0, early-type spirals, late-type spirals, and 
irregulars). In addition, we want to compare the properties 
of galaxies in the Virgo cluster to those in the field. We esti- 
mated that 18 galaxies per bin were required to obtain good 
statistics on the average properties of each bin, or 144 galax- 
ies total. However, since our sample only contains 148 Virgo 
galaxies, we compromised on 9 galaxies per bin in Virgo and 
18 per bin in the field, giving a sample size of 108 galaxies. 
In addition, we limited the field and Virgo samples to galax- 
ies with D25 < 4' so that we could use the jiggle-map mode 
with HARP-B (|Buckle et all [2009]) and still map the inner 
quarter of each galaxy. 

To select the final sample, the field and Virgo lists 
were each divided into the four morphological bins de- 
scribed above using the numerical Hubble stage, T (E/S0: 
T ^ 0.5; early-type spirals: 0.5 < T ^ 4.5; late-type spirals: 
4.5 < T < 9.5; irregulars: T > 9.5). For the field samples, 
we first applied an H I flux cutoff of 6.3 Jy km s _1 and then 
randomly selected 18 galaxies in each of the four morpholog- 
ical bins (for the elliptical galaxies, there were precisely 18 
galaxies above this cutoff). The smaller Virgo sample meant 
that we could only use this same process for the late-type 
spiral bin. The brightest 9 in each category of E/S0 (^ 2.1 
Jy km s _1 ), Irr (^ 4.1 Jy km s _1 ), and early-type spiral 
galaxies 6.2 Jy km s _1 ) from Vir go were included, in ad- 
dition to a random selection of 9 late-type spirals with H I 
fluxes brighter than 6.3 Jy km s . The final field and Virgo 
source lists are given in Tables IATTA31 in Appendix [£] 

2.2 CO J=3-2 observations and data processing 

CO J=3-2 observations for all galaxies were obtained on 
the JCMT between 2007 November and 2009 November. 
The angular resolution of the JCMT at this frequency is 
14.5", which corresponds to a linear resolution ranging from 
0.2 to 1.2 kpc for the galaxies in our sample. The 30 large 
(D25 > 4') SINGS galaxies in the NGLS were observed in 
raster map mode, w hile the smaller gal axies were observed 
in jiggle-map mode dBuckle et al.l I2009T ) . All galaxies were 
mapped over a rectangular area corresponding to D25/2 on 
a side (with raster maps oriented with the position angle 
of the galaxy's semi-major axis), with a 1 sigma sensitivity 
of better than 19 mK (T\) (32 mK T M b for t)mb = 0.6) 
at a spectral resolution of 2 km s~ . We use d the 16 
pixel array receiver HARP-B (|Buckle et al.|[2009l) with the 
Auto-Correlation Spectral Imaging System (ACSIS) corre- 
lator configured to have a bandwidth of 1 GHz and a resolu- 
tion of 0.488 MHz (0.43 km s" 1 at the frequency of the CO 
J=3-2 transition) . The total time allocated to the CO J=3- 
2 portion of the NGLS, including pointing and calibration 
observations, was 256 hours. 

Calibration was checked each night by observing one or 
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Table 1. Large (D25 > 4') galaxies from the SINGS sample 
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+69:40:47 


5.2 


2.2 


65 


254 


NGC 3184 d 


10:18:17.0 


+41:25:28 


4.0 


3.8 


135 


574 


NGC 3198 d 


10:19:55.0 


+45:32:59 


3.9 


1.3 


35 


666 


IC 2574 d 


10:28:23.5 


+68:24:44 


6.7 


3.2 


50 


28 


NGC 3351 d 


10:43:57.8 


+11:42:13 


3.8 


2.2 


13 


768 


NGC 3521"* 


11:05:48.9 


-00:02:06 


5.4 


2.7 


163 


778 


NGC 3627 d 


11:20:14.9 


+12:59:30 


1.4 


2.0 


173 


697 


NGC 3938 d 


11:52:49.3 


+44:07:17 


2.5 


2.4 





829 


NGC 4236 


12:16:42.0 


+69:27:46 


11.4 


3.6 


162 


10 


NGC 4254 d 


12:18:49.6 


+14:25:00 


2.7 


2.3 





2412 


NGC 4321 d 


12:22:54.9 


+15:49:21 


3.8 


3.1 


30 


1599 


NGC 4450 


12:28:29.7 


+17:05:06 


2.6 


1.9 


175 


1905 


NGC 4559 


12:35:57.7 


+27:57:36 


6.0 


2.1 


150 


827 


NGC 4569 d 


12:36:50.0 


+13:09:46 


5.2 


2.3 


23 


-179 


NGC 4579 d 


12:37:44.0 


+11:49:07 


2.8 


2.2 


95 


1540 


NGC 4594 


12:39:59.3 


-11:37:22 


4.2 


2.0 


90 


1088 


NGC 4631 d 


12:42:07.7 


+32:32:34 


7.4 


1.3 


86 


640 


NGC 4725 


12:50:26.7 


+25:30:03 


5.5 


3.7 


35 


1188 


NGC 4736 


12:50:53.1 


+41:07:12 


6.3 


5.5 


105 


295 


NGC 4826 


12:56:43.9 


+21:41:00 


5.1 


2.6 


115 


390 


NGC 5033 d 


13:13:27.6 


+36:35:38 


5.1 


1.9 


170 


873 


NGC 5055 d 


13:15:49.4 


+42:01:46 


6.6 


4.0 


105 


495 


NGC 5194 d 


13:29:52.4 


+47:11:41 


4.9 


3.4 


163 


452 



a Prom RC3 (de Vaucouleurs et al. 1991) 
b From Buta et al. (2007) 

c Systemic velocity from the H I line (heliocentric) 

d Member of high-priority SINGS subset for SCUBA-2 observations 



more of a set of six calibrator sources in the CO J=3-2 
transition. The calibration sources used were W75N, CRL 
618, IRC+10216, NGC 2071IR, CRL 2688, and IRAS 17293- 
2422. The peak intensities of the calibrator spectra agreed 
to within ±8%, while the integrated intensities agreed to 
within ±13% from 2007 November to 2008 July and to ±9% 
from 2008 August onwards. We adopt 10% as the internal 
calibration uncertainty on the NGLS. Pointing was checked 
before starting a new source, and every 1-2 hours (more 
frequently near sunrise and sunset). Any scan where the 
pointing offset that was measured after the scan differed by 
more than 4" from the pointing measured at the start of the 
scan was rejected in the final data analysis. Less than 5% 
of the data were rejected because of large pointing changes. 
The rms of the pointing measurements on a given night was 
typically better than 2". 

Two sources in the NGLS, NGC 5055 and NGC 337, 
were observed during Science Verification runs in 2007 May 
and 2007 July. For these observing runs, the strength of the 
spectral line calibrators was only 85% of the expected values 
and so the data for these two galaxies have been rescaled by 
a factor of 1.18 to bring them onto the same calibration scale 



as the rest of the survey. NGC 4789A (DDO 154) was ob- 
served as part of a separate program (M07AC14, PI B. War- 
ren), as was NGC 5194 (M51) (M06AN05, PI R. Tilanus). 
For both these galaxies, the integration time per point in 
the map is longer than that of a typical NGLS observation. 

Details of the r eduction of the raster map data are given 
in the appendix in IWarren et alJ l)2010t ) and so we discuss 
only the jiggle map processing in detail here. The individ- 
ual raw data files were flagged to remove data from any of 
the 16 individual receptors with bad baselines and then the 
scans were combined into a data cube using the "nearest" 
weighting function to determine the contribution of individ- 
ual receptors to each pixel in the final map. This weighting 
function causes each measurement from a raw data file to 
contribute to only a single pixel in the final cube. The pixel 
size in the maps is 7.5". We fit and subtracted a first-order 
baseline from the data cube while excluding the central 400 
km s _1 region from the baseline fit. 

We used the clumpfind algorithm dWilliams et al.|[l994l ) 
implemented as part of the Clump Identification and Anal- 
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Table 2. Small (D 25 < 4')galaxies from the SINGS sample 



IN EL III C 


T9nnn f\\ a 


R( T9non c\\ a 


Vu , b 
yhel 


AiXGrilELtlVC IN EilllG 




(h.m.s) 


( ■ ■ ) 


(km s ) 




NGC 0337 


00:59:50.1 


-07:34:41 


1664 




NGC 0584 


01:31:20.8 


-06:52:05 


1861 




NGC 0855 


02:14:03.6 


+27:52:38 


583 




PGC 023521 


08:23:55.0 


+71:01:57 


113 


M81DwA 


UGC 05139 


09:40:32.3 


+71:10:56 


137 


Hoi 


NGC 3049 


09:54:49.7 


+09:16:18 


1458 




UGC 05336 


09:57:32.1 


+69:02:46 


48 


HoIX 


UGC 05423 


10:05:30.4 


+70:21:54 


344 


M81DwB 


NGC 3190 


10:18:05.3 


+21:49:58 


1306 




NGC 3265 


10:31:06.8 


+28:47:47 


1410 




UGC 05720 


10:32:32.0 


+54:24:04 


1456 


Mrk33 


NGC 3773 


11:38:13.0 


+12:06:43 


1002 




NGC 4625 


12:41:52.7 


+41:16:26 


590 




NGC 4789A 


12:54:05.5 


+27:08:55 


374 


DD0154 


UGC 08201 


13:06:24.9 


+67:42:25 


34 


DDQ165 


NGC 5474 


14:05:01.6 


+53:39:45 


240 





a Prom RC3 fde Vaucouleurs et al. 1991) 

b Systemic velocity from the H I line (heliocentric) 



ysis Package (CUPIE0) (|Berrv et alll2007h task findclumps 
to identify regions with emission above a specified signal- 
to-noise ratio. The algorithm measures the mean noise in 
the data cube as part of the processing step. However, the 
sensitivity of each of the sixteen individual receptors in the 
HARP-B array varies quite significantly, fn addition, a por- 
tion of the data were obtained when four of the receptors 
were inoperative and a rotation of the K-mirror between 
two observations was used to obtain a complete map of the 
desired area. As a result, the noise in a final data cube can 
vary quite significantly from region to region in a jiggle map. 
Thus, the mean noise that findclumps measures will under 
or over-estimate the noise in different regions of the cube. 

To mitigate the effect of this noise variation, we calcu- 
lated an image of the noise in the map by averaging the noise 
in the line free regions at the end of the spectrum. We then 
divided the cube by this noise image to produce a cube which 
was now in units of signal-to-noise, rather than Kelvins. This 
data cube was then boxcar smoothed by 3 pixels and 25 ve- 
locity channels. We then applied the findclumps algorithm 
to this signal-to-noise cube to identify regions with emis- 
sion with signal-to-noise greater than 2a, 2.5a, or 3a. In all 
cases, the spacing between contour levels (the AT paeram- 
eter) was set to 2a. The mask of regions of real emission 
produced by findclumps was then applied to the original, 
unsmoothed data cube and this masked cube was then col- 
lapsed in velocity to produce moment maps. For the raster 
maps, we produced masks with signal-to-noise cutoffs of 2a, 
2.5a, and 3a. 

The zeroth moment map measures the integrated in- 
tensity in the data cube, J Tdv = ^TiAv where Av is 
the channel width, Ti is the temperature in an individual 
channel. To obtain the final integrated intensity map, the 

c CUPID is part of the Starlink iCurrie et alj 12008ft 
software package, which is available for download from 
http: / /starlink.jach. hawaii.edu 

© 0000 RAS, MNRAS 000, 000-000 



zeroth moment map needed to be multiplied by the noise 
map (to recover units of K km s _1 ) and also to be divided 
by t\mb = 0.6 to convert to the Tmb temperature scale. 
For the zeroth moment maps, we use a noise cutoff of 2a; 
we chose to use this cutoff after comparing the CO fluxes 
obtained with different noise cutoffs (Figure [TJ . For bright 
galaxies with a flux greater than 150 K (T^) km s _1 pixels, 
the data show that fluxes measured using a 2a cutoff are sys- 
tematically offset from fluxes measured with a 3a cutoff, in 
the sense that the maps made withe the 3a cut off underesti- 
mate t he total flux. Tw o galaxies, NGC 2403 (Bendo et al.l 



l2010ah . and NGC 3031 (|Sanchez-Gallego et al][201ll h have 



much larger offsets and were not included in deriving the 
average offset. These galaxies both have such low CO sur- 
face brightnesses such that an increase in the noise cutoff 
causes large areas of the galaxy to be masked out and thus 
has a large impact on the total flux that is measured. 

The first moment map measures the mean velocity of 
the emission v = ^ TiVi/^2 Ti where v% is the velocity of a 
given velocity channel. The mean velocity calculated in this 
way can be rather sensitive to noise spikes, and so we used 
a higher signal-to-noise threshold of 2.5a in calculating the 
velocity field maps. 

The second moment maps measure the velocity disper- 
sion, a v , for each pixel in the image using 

a v = N^ Ti(vi — v) 2 / N^ Ti (1) 

This method of calculating the velocity dispersion differs 
fr om those used in prev i ous ex tragalactic CO studies such 
as lCombes fc Becauaertl (|l997T ). which typically fit gaussian 
profiles to the CO lines. However, in the limit of gaussian 
lines with a high signal-to-noise ratio, the values from the 
moment 2 maps should agree with the results from fitting a 
gaussian directly to the line profiles. Note that for a gaussian 
line profile, a v is equal to the full-width at half maximum 
of the line divided by a factor of 2.355. The second moment 
maps were calculated using a signal-to-noise threshold of 
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2000 




500 1000 1500 2000 

CO J = 3-2 flux, 2a cutoff (K km s _1 pixels, T s *) 

Figure 1. CO J=3-2 flux for galaxies in the NGLS compared 
for a 2a cutoff versus a 3a cutoff in making the moment maps. 
The solid line represents 2a flux = 152 + 3a flux (see text). NGC 
3034 (M82) is much brighter than all other galaxies and is not 
included in this plot; its flux offset is consistent with the rest of 
the sample. Galaxies to the left of the vertical dashed line were 
not included in calculating the average offset (see text). 

2.5a. A more detailed discussion of the possible systema tics 
in second moment maps is given in IWilson et all (|201ll ). 

CO J=3-2 images for the SINGS sample are shown in 
Appendix [C] Images for the Virgo and field samples will be 
presented in a future paper (Golding et al., in prep.). The 
reduced images, noise maps, and spectral cubes are available 
via the survey websit^fl and will eventually be made available 
via the Canadian Astronomical Data Centre (CADC). We 
note that some of the images show residual noise effects, par- 
ticularly around the outer edges of the map (see NGC 4736 
and NGC 5033 for particularly obvious examples). These ar- 
tifacts are generally due to the increasing noise level in the 
data cubes towards the edges of the map, which due to the 
scanning technique do not have as long an integration time 
and are not sampled by as many of the individual HARP 
detectors. Anyone interested in the reality of a particular 
faint feature is encouraged to consult the images and noise 
maps available on the survey web site. 

2.3 Future SCUBA-2 observations 

The Submillimetre Common-User Bolometer Array-2 
(SCUBA-2, Holland et al. 2006) is a large-format bolometer 
camera for the JCMT that is designed to produce simulta- 
neous continuum images at 450 and 850 /im. The camera 
consists of four 32 x 24 transition-edge sensor arrays for 
each of the two wavelengths and has a total field of view 
of ~ 7' x 7'. The SCUBA-2 portion of the NGLS began 
in December 2011 with a period of science verification ob- 
servations for the legacy surveys. In this section, we briefly 
describe our planned observing strategy with SCUBA-2. 

The ultimate goal of the survey is to map all the galax- 
ies in the NGLS with SCUBA-2 out to D 25 to a la limit of 
1.7 mjy at 850 /im; this sensitivity limit includes both in- 
strumental noise and confusion noise from the high-redshift 

d http:/ /www. physics. mcmaster.ca/~wilson/www_xfer/NGLS/ 



extragalactic background. Good 450 /im data provide im- 
proved spatial resolution and are critical to tracing the up- 
turn of the dust spectral energy distribution that indicates 
the presence of very cold dust or a change in dust properties 
(see § [T]). With similar (simultaneous) integration times at 
both 450 and 850 /j,m, we estimate that the 450 /im sensi- 
tivity will be about a factor of two worse than at 850 fim for 
a given surface density of dust and smoothed to the same 
angular resolution. 

To put these sensitivities in physical terms, we can con- 
vert to an equivalent mass surface density (gas plus dust) by 
using the formula given in Johnstone et al. (2001). Assuming 
Tdust = 20 K, K850 = 0.0043 cm 2 g _1 , and a circular beam 
with diameter 15", a 4a detection at 850 /im corresponds 
to an average visual extinction of A v — 3 mag, which corre- 
sponds to surface densities of 6 x 10 21 H/cm 2 or 45 M /pc 2 . 
For a galaxy at a distance of 10 Mpc, the equivalent mass 
sensitivity is 2 x 10 7 Mg, equal to a few massive giant molec- 
ular clouds within a 730 pc diameter region. 



3 CO PROPERTIES 

3.1 Measuring CO luminosities 

For galaxies which were clearly detected, the CO J=3-2 lu- 
minosity was measured directly from the moment images 
made with a 2a cutoff mask. We used apertures chosen by 
eye to capture all the emission from the galaxy while ex- 
cluding the occasional noisy or doubtful emission towards 
the edges of the map. Thus, the CO luminosity is given by 

Lco(3-2) = ^2 H 00 ^ x (4.848D Mpc X L ptx f (2) 

where I(CO)i is the CO intensity in an individual pixel of 
the map, Dm P c is the distance to the galaxy in Mpc, and 
Lpix is the size of an individual pixel in the map (7.5" for 
jiggle maps and 7.2761" for raster maps). 

We calculated a map of the uncertainty in I (CO) for 
each pixel, a;, as 

01 = &chanV Av c hanAvi ine y / 'l + Avu„e / Avbase (3) 

where Av c han is the velocity width of a single channel, a c fcan 
is the standard deviation in K of the line-free channels, 
Avune is the velocity range used to measure the line, and 
Avbase is the velocity range used to fit the baseline. Now, 
both the raster and the jiggle observing modes use a shared 
off position observed for a longer integration time than each 
of the individual on positions. Because the off position is 
shared across many pixels in the map, a simple combination 
of the noise for each pixel in the map will underestimate the 
true uncertainty. Assuming an on-source integration time 
per pixel of t on and a corresponding integration time on the 
off position of t ff, the uncertainty in the CO luminosity is 
given by 

a L{CO ) = \J(1 - a) ^ a l + a C^2 cr i) 2 x( 4 -848D M pcXLp^) 2 

where a = l/(t ff/t n + 1) (see Appendix |B)) . For the jiggle 
map mode, t ff = 4i on and so a = 0.2. For the raster maps, 
the ratio of t f / /t on varied with the size and dimensions of 
the map; we have estimated an average value a — 1.07 for 
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Table 3. CO and far-infrared luminosities for the SINGS sample 



Name 


D a 


AT 6 


T c 

^CO(3—2) 


lot? TjTTT t> 


^ r 1 til "00(3' 


-2) 






(Mac) 
V 1VA PW 


(mK) 


(10 7 K km s _1 pc 2 ) 


(Ln) 


L Q /(K km s _1 


PC 2 ) 




NGC 0024 


7.6 


27 


< 0.4 


8.19 


> 40 




< 0.10 


NGC 0628 


7.3 


24 


5.2 ± 1.0 


9.55 


67 




0.25 


NGC 0925 


9.1 


25 


0.9 ± 0.2 


9.31 


221 




0.04 


NGC 2403 


3.2 


23 


1.7 ± 0.3 


9.09 


72 




0.12 


UGC 04305 


3.4 


31 


< 0.09 


7.51 


> 38 




< 0.02 


NGC 2841 


14.1 


24 


< 1.1 


9.45 


> 249 




< 0.02 


NGC 2976 


3.6 


26 


0.50 ± 0.08 


8.72 


103 




0.65 


NGC 3031 


3.6 


26 


0.91 ± 0.35 


9.41 


284 




0.05 


NGC 3034 


3.6 


31 


39.2 ± 0.4 


10.60 


102 




8.6 


NGC 3184 


11.1 


28 


9.5 ± 1.4 


9.61 


43 




0.55 


NGC 3198 


13.7 


31 


6.9 ± 1.1 


9.61 


60 




0.12 


IC 2574 


3.8 


23 


< 0.08 








< 0.01 


NGC 3351 


9.3 


23 


5.1 ± 0.6 


9.66 


88 




0.89 


NGC 3521 


7.9 


27 


17.9 ± 1.2 


9.97 


52 




0.73 


NGC 3627 


9.4 


23 


31.1 ± 1.7 


10.18 


19 




6.5 


NGC 3938 


14.7 


26 


12.3 ± 2.0 


9.85 


57 




0.51 


NGC 4236 


4.4 


21 


0.53 ± 0.15 


8.19 


29 




0.05 


NGC 4254 


16.7 


27 


73.8 ± 4.8 


10.50 


43 




2.7 


NGC 4321 


16.7 


28 


54.7 ± 5.3 


10.36 


42 




3.1 


NGC 4450 


16.7 


23 


1.1 ± 0.3 


9.06 


106 




0.62 


NGC 4559 


9.3 


23 


1.9 ± 0.4 


9.42 


138 




0.07 


NGC 4569 


16.7 


26 


20.1 ± 1.7 


9.95 


14 




5.3 


NGC 4579 


16.7 


20 


7.9 ± 1.4 


9.84 


87 




2.3 


NGC 4594 


9.8 


28 


< 0.7 


8.92 


> 122 




< 0.43 


NGC 4631 


7.7 


20 


14.6 ± 0.7 


10.10 


85 




0.31 


NGC 4725 


11.9 


25 


< 0.9 








< 0.05 


NGC 4736 


5.2 


30 


5.2 ± 0.4 


9.65 


86 




1.9 


NGC 4826 


7.5 


30 


9.7 ± 0.5 


9.75 


58 




3.3 


NGC 5033 


16.2 


24 


23.1 ± 3.0 


10.19 


67 




1.6 


NGC 5055 


7.9 


30 


19.7 ± 1.7 


10.00 


51 




0.63 


NGC 5194 


7.7 


15 


53.2 ± 2.0 


10.21 


31 




4.0 


NGC 0337 


23.1 


37 


< 5.6 


10.13 


> 240 




< 0.15 


NGC 0584 


20.1 


19 


< 2.5 








< 2.6 


NGC 0855 


9.7 


15 


< 0.5 


8.29 


> 41 




< 0.54 


PGC 023521 


3.6 


29 


< 0.1 








< 1.4 


UGC 05139 


3.6 


22 


< 0.09 








< 0.12 


INCjC 3049 


22.7 


20 


1.4 i 0.3 


9.39 


176 




0.15 


UGC 05336 


3.6 


21 


< 0.09 








< 0.0.05 


UGC 05423 


3.6 


24 


< 0.1 








< 1.5 


NGC 3190 


21.3 


19 


< 2.7 


9.52 


> 121 




< 0.86 


NGC 3265 


22.7 


20 


< 3.3 








< 2.7 


UGC 05720 


25.0 


21 


2.8 ± 0.6 








0.70 


NGC 3773 


10.5 


16 


0.20 ± 0.05 


8.38 


122 




0.54 


NGC 4625 


9.4 


19 


0.12 ± 0.03 


8.37 


196 




0.03 


NGC 4789A 


3.8 


7.3 


< 0.03 








< 0.02 


UGC 08201 


4.6 


18 


< 0.1 








< 0.19 


NGC 5474 


6.5 


26 


< 0.3 


8.14 


> 40 




< 0.06 



a NGC 925, NGC 2403, NGC 3031, NGC 3198, NGC 3351, NGC 3627, NGC 4725: Freedman et al. (20011 ; NGC 2976, NGC 4236, 

UGC04305, UGC08201: Karachentsev et al. (2002) NGC 3034, PGC 023521, UGC05139, UGC05336, UGC 05423: same distance 

as NGC 3031 fFreedman et al. 2001) : NGC 584, NGC 855, NGC 4594, NGC 4736, NGC 4826, NGC 5194: Tonrv et al. (2001) ; 

Virgo Cluster: Mei et al. (2007) : NGC 628: Karachentsev et al. (2004) : IC 2574: Dalcanton et al. (2009) NGC 2841: 

Macri et al. (2001) : NGC 3184: Leonard et al. (2002) : NGC 4631: Seth et al. (2005) . Remaining galaxies: Hubble flow distance with 

velocity corrected for Virgo infall (Mould et al. 2000) and H a = 70.5 km s — 1 Mpc -1 . 

b rms noise in individual spectra in the data cube at 20 km s 1 resolution on Tmb scale. 

c Upper limits are 2<r limits calculated over an area of 1' and a line width of 100 km s — . 

d Lpjn from^anders_et_aL (_2003) and ljis^nfekd^t_ah Q2007) (see text) adjusted for distances given here. 

e Total Hi masses (including emission from areas outside the CO disk) from fluxes in Walter et al. (2008), 

Chung^t_aL (2009), or the HyperLeda database (Paturel et al. 2003); see text. Mjj 2 calculated assuming a CO J=3-2/J=l-0 

line ratio of 0.18, which may overestimate the total mass in some galaxies, as well as a standard value for Xco\ see text. 
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the rasters and have used this in the noise calculations for 
all the galaxies. 

For galaxies without detections, we calculated an uncer- 
tainty map using equation (3) by adopting a line width of 
100 km s _1 and then calculated ol(CO) using equation (4). 
For galaxies observed in jiggle map mode, the velocity width 
for the baseline determination was always 400 km s _1 . For 
the galaxies observed in raster mode, the actual region used 
to determine the baseline varied slightly from pixel to pixel 
as a result of weak emission appearing from place to place in 
the spectrum. For these galaxies, we used the actual region 
used for the baseline in calculating <Xj for each pixel. 

Of the 31 large galaxies in the SINGS sample, 25 galax- 
ies were detected for a detection rate of 77%. Of the 16 
smaller galaxies, only four galaxies (NGC3049, NGC 3773, 
NGC 4625, and UGC 05720) are detected at > 4<r in the 
integrated intensity maps. A fifth galaxy, NGC 3190, is 
marginally detected at the 3cr level. Thus, the detection rate 
for the small SINGS galaxies is 25%, significantly lower than 
for the larger SINGS galaxies. These smaller SINGS galaxies 
tend to be physically smaller galaxies (often dwarf galaxies) 
and also have on average lower Hi fluxes than their larger 
counterparts, both of which may impact the detection rate. 



3.2 Comparison with other surveys and the global 
CO line ratio 

Because the different CO lines have different excitation tem- 
peratures and critical densities (JEJ, comparing the emission 
from different CO lines can provide information on the phys- 
ical conditions in the molecular gas. These different excita- 
tion conditions may also produce different radial distribu- 
tions in the different CO lines . The HERA CO- Line Extra- 
galactic Survey, HERACLES, l|Lerov et alj2009l ) mapped 18 
nearby galaxies in the CO J=2-l line. The large area cov- 
ered per galaxy and excellent sensitivity make the data from 
this survey an excellent comparison to the NGLS. Figure [2] 
compares the CO luminosities for the 9 galaxies that were 
detected in both surveys. (Three additional galaxies in com- 
mon had only upper limits in both surveys, and NGC 2841 
was detected by HERACLES but not by us.) The correlation 
between the two sets of global measurements is good, with 
a mean CO J=3-2/2-l line ratio of 0.36 ± 0.04 (where the 
quoted uncertainty is the stan dard deviation of t he mean; 
the standard deviation is 0.13). iKuno et all (|2007l ) mapped 
40 nearby galaxies in the CO J=l-0 line, also covering a large 
area with good sensitivity, although the fact that the maps 
are somewhat undersampled may introduce some scatter. 
Figure [2] compares the CO luminosities for the 11 galaxies 
that were detected in both surveys. The correlation between 
the two sets of measurements is again good, with a mean CO 
J=3-2/l-0 line ratio of 0.18±0.02 (standard deviation 0.06). 

This average CO J=3-2/l-0 line ratio is somewhat 
smaller than the line ratio s (0.4-0.8) seen in i ndividual giant 
molecular clouds in M33 I Wilson et al.lll997fl and is at the 
low end of th e range (0.2-1.9) mea s ured in the central regions 
of galaxies (iMauersbergerl 1 19991; iMao et al.l 120101 ). How- 
ever, the sample of Mao et alT i 201oh includes roughly 50% 
Seyfert, LINER, merging, or OH megamaser galaxies which 
tend to show a higher than average line ratio. Not surpris- 
ingly, our average line ratio is also at the small end of the line 
ratios (0.1-1.9) measured for low-redshift luminous infrared 



galaxies (|Leech et al.ll201ol ; |Papadopoulos et al.ll2~012T ). Note 
also that we have not attempted in either comparison to 
correct for differences in the fraction of the galaxy mapped 
in each survey. Given the differences in the processing tech- 
niques adopted between the three surveys, the small scatter 
in this relationship is a good indicator that both surveys are 
successfully measuring global CO luminosities in the low 
signal-to-noise regime. A more detailed analysis of spatially 
resolved images of the CO line ratios by combining these 
various surveys with new 13 CO data will be presented in 
Rosolowsky et al. (in prep.). 



4 COMPARISON WITH OTHER GLOBAL 
PROPERTIES 

4.1 The molecular mass fraction in the 
interstellar medium 

To examine the molecular gas mass fraction in the inter- 
stellar medium of the SINGS galaxies in our sample, we 
have converted the CO J=3-2 luminosities to molecular 
hydrogen mass adopting a CO-to-H2 conversion factor o f 
X CO = 2 x 10 20 cm" 2 (K km s -1 ) -1 (|Strong et alJll988h . 
We have adopted a CO J=3-2/J=l-0 line ratio of 0.18 
whi ch is the mean va lue derived from our comparison with 
the lKuno~e t al. (2007) sample. With these assumptions, the 
molecular hydrogen mass is given by 



H 2 



17.8(^/0.18)" 



J CO(3-2) 



(•») 



where f?3i is the CO J=3-2/J=l-0 line ratio, Mh 2 is in M©, 
and Lco{3-2) is in units of K km s _1 pc 2 . We note that 
the mean value for the line ratio may not be appropriate 
for all galaxies; in particular, it may result in an overes- 
timate of the H2 mass for more luminous galaxies and an 
underestimate in less luminous galaxies. We also note that 
adopting a single value for the CO-to-H2 conversion factor, 
Xco , will likely result in an underestimate of the molecular 
gas mass in galaxies or regions of galaxies where th e metal- 
licity is more than about a factor of two below solar (|Wilsonl 
ll995l ; lArimoto et al.|[l99r3 ; lBolatto et al.ll2008h . Finally, our 
adopted value of Xco will likely result in an overestimate 
of the gas mass in the starburst galaxy M82, as lower val- 
ues of Xco are m ore appropriate in starbur st and luminous 
infrared galaxies (|Downes fc SolomorJ[l99sl ). 

We have compiled H I fluxes from the lite rature for all 
of our samp le. The prefe re d ref erences were IWalter et al.l 
Hooi) and IChung et all j2009h . with most of the re- 
maining galaxies re trieved from the HyperLeda database 
(iPaturel et all l2003h . The Hi flux for UGC 08201 is from 



ICannon et al 



(120111). wh ile the Hi flux for NGC 3190 is 
taken from lMartinl l|l998). All fluxes have been converted to 
Hi mass es using the d is tance s in Table [3] and the equations 
given in IWalter et all (|200Sl V Figure [3] shows the H2/H1 
mass fraction as a function of Hi mass for all the SINGS 
galaxies. Of the galaxies with CO detections, 10 galaxies 
have H2/H1 mass fractions greater than 1, 9 galaxies of 
0.5 < Mh 2 /Mhi < 1 and 10 galaxies have H2/H1 mass frac- 
tions less than 0.5. The 18 galaxies for which we have only 
CO upper limits span a similar range of atomic gas masses, 
although 78% of these galaxies have Hi masses less than 
1.5 x 10 Mq. On average, these galaxies have significantly 
lower molecular gas fractions than the detected galaxies. 
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Figure 2. (top) LcO(3-2) versus ^co{2— 1) f° r all galaxies detected in both the NGLS and the HERACLES survey llLerov et al.ll2009h . 
The HERACLES luminositi es have been adjus ted to use the distances adopted in this paper. Calibration uncertainties in the HERACLES 
data are estimated at 20% jLerov et al.ll2009l ). The line shows the mean CO J=3-2/2-l ratio of 0.36 determi ned from the data points 
(see text), (bottom) LcO(3-2) versus L<;jo(i-o) for all galaxies detected in both the NGLS and the survey of lKuno et al ] d2007ft . The 
CO J=l-0 luminosities have been adjusted to use the distances adopted in this paper. The line shows the mean CO J=3-2/l-0 ratio of 
0.18 determined from the data points (see text). 



o 



q n 



The H2/H1 mass fractions shown here are on average 
somewhat l arger than s i milar measurements for 14 galax- 
ies given in iLerov et al.l (|2009l ). For the 9 gala xies in com- 
mon w ith our survey, the mass fraction given in lLerov et al.l 
(2009) is a factor of ~ 1.5 times smaller than the values for 
the same galaxies given in Table [3] The primar y reason for 
this d ifference appears to be the CO line ratio. ILerov et al.l 
(2009) adopt a CO J=2-l/J=l-0 line ratio of 0.8 measured 
from comparing the pea k line temper a tures with several sur- 
veys, including that of iKuno et all (|2007T> . On the other 
hand, comparing the global CO luminosities for five galax- 
ies in common between these two surveys gives an average 
CO J=2-l/J=l-0 ratio of 0.52 ±0.06 . If we rescale the mass 
fractions given in ILerov et all |2009h to this lower line ratio, 
the average agreement between their mass fractions and our 
values is good. 

This comparison illustrates the importance of using the 
appropriate value of the CO line ratio when observing galax- 
ies in the higher CO transitions. It is likely too simplis- 
tic to calculate gas masses using a single value of the line 
ratio when the galaxy properties vary widely. Indeed, the 
molecular gas mass calculated for NGC 3034 (M82) is cer- 
tainly an overestimate because the CO line ratios in this 



starburst galaxy are much larger (Ward et al. 20031 ) . Simi- 
larly, at least two of the bright spirals in Virg o (NGC 4321 
and N GC 4254) also show elevated line ratios (|Wilson et all 
2009). In contrast, for low luminosity galaxies the molecu- 
lar gas mass may be underestimated if the CO excitation 
is very low. NGC 2841 ha s been detected in the CO J=2-l 
line by ILerov et~aH l|2009h ; adopting a CO J=2-l/l-0 line 
ratio of 0.52 leads to a molecular mass fraction of 0.20, a 
factor of 10 above the upper limit deduced from the CO 
J=3-2 upper limit. The fact that we have measured upper 
limits to the CO luminosity over a 1' diameter aperture, 
which is roughly a factor of 7 smaller than the area of the 
galaxy within D25/2, may go some way to explaining this 
discrepancy. 

4.2 The Lqo(3-2)—Lfir correlation 

We have used the NG1S CO measurements to examine the 
correlation betwe en the far-infrared luminosity and the CO 
J=3-2 luminosity. Ilono et all (|2009l ) presented a nearly lin- 
ear correlation between these two luminosities, and sug- 
gested that the CO J=3-2 line could be a good tracer of the 
dense gas involved in star formation, similar to the HCN 
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Figure 3. The ratio of H2 to Hi mass plotted as a function of H I mass. Galaxies with only upper limits to the CO luminosity arc shown 
as arrows. H2 masses have been calculated assuming a CO J=3-2/J=l-0 line ratio of 0.18; this line ratio may not be appropriate for all 
galaxies, particularly the more luminous ones such as NGC 3034 (M82). 
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Figure 4. Lpjn plotted versus l^cO('A-2) f° r an SINGS galaxies in the NGLS with far-infrared data. Galaxies which were not detected 
in CO are indicated by arrows. 
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Figure 5. Lfir/ Lqq^— 2) plotted versus ^co{3—2) f° r an SINGS galaxies in the NGLS with far-infrared data. Galaxies which were 
not detected in CO are indicated by arrows. 



J=l-0 line (|Gao fc Solomon! 12004 ). We have collected far- 
infrared luminosities derived from IRAS data for all the 
SINGS galaxies for which such measurements are available. 
The preferred sourc e is the Rev ised Bright Galaxy Su rvey 
l|Sanders et alj|2003ft followed bv lLisenfeld et alj i|2007l ). For 
8 additional galaxies we retrieved the 60 and 100 fim fluxes 
from the IRAS Faint Source Catalog and calculated the 
far-inf rared luminosity using the formula in lLisenfeld et ail 
(|2007t ). When necessary, published far-infrared luminosities 
were adjusted for our adopted distances. 

Figure [4] plots the far-infrared luminosity as a function 
of the CO J=3-2 luminosity, while Figure [S] plots the ratio 
Lfir/ Lco(3~2) as a function of the CO luminosity. These 
figures reveal a remarkably tight correlation for the more 
luminous galaxies and increased scatter among the fainter 
galaxies. Of 19 galaxies with 9.5 < log Lfir < 10.7, only 
NGC 337 is not detected in the CO J=3-2 line. In com- 
parison to NGC 2841, which has a similar lower limit to 
Lfi r/Lcq(3-2) and w hich has been detected in the CO 2-1 
line (|Lerov et al.ll2009l ). NGC 337 is five times more lumi- 
nous in the infrared yet with a somewhat smaller atomic gas 
mass and linear extent (as measured by th e D25 diameter). 
Desp ite its inclusion in the SINGS sample (|Kennicutt et al] 
l2003t ). this galaxy has not been very well studied. One ex- 
planation of its unusually high infrared luminosity would 
be t he presence of an AG N. However, the GALEX im- 
ages l|Gil de Paz et all 1201 ll ) show no hint of a bright nu- 
clear source while the mid-infra red emission is extended and 
asymmetric (|Bendo et al.ll2007h . Its optical through infrared 
spectral en ergy distribution is consistent with a normal spi- 
ral galaxy l|Dale et al.l I2007T ) . optical spectroscopy is dom- 
inated by star formation indicators, and the metallicity is 



near solar (|Moustakas et al1 l2010). This galaxy remains a 
bit of a mystery and would be worthy of further study. 

Excluding NGC 337, the mean Lfir/Lqo(3-2) ratio 
for the luminous galaxies galaxies is 62 ± 5 with a standard 
deviation of 20. The situation is more complicated for the 
fainter galaxies. Of 12 galaxies with 8.3 < log Lfir < 9.5, 
only 7 (58%) have global CO luminosities that are detected 
at better than the 4<r level. For these 7 galaxies, the mean 
Lfir/ Lco(3-2) is 147 ± 20 with a standard deviation of 
53. Some of the increased scatter is likely due to the larger 
CO measurement uncertainties on these fainter galaxies, but 
this seems unlikely to explain all the scatter and the in- 
creased mean ratio. One possibility is that these lower lu- 
minosity and hence likely lower mass galaxies may tend to 
have lower metallicities, in which case the CO luminosity 
may systematically underestimate the molecular hydrogen 
gas mass and cause them to ap pear underluminous in CO 
(israelii 19971 ; iLerov et al.l [201ll ). However, the characteris- 
tic oxygen abundance for these lower luminosity galaxies is 
not signi ficantly different from t hat of the higher luminosity 
systems IjMoustakas et alJl201CT ). Another possibility is that 
these fainter galaxies have lower average CO surface bright- 
nesses, in which case we could be systematically underes- 
timating the CO luminosity due to the low signal-to-noise 
ratio in the data. An example of this type of effect is seen for 
NGC 2403, for which the CO luminosity increases by more 
than a factor of 2 when the moment maps are made with 
a 2a cutoff compared to a 3a cutoff (see Section 12. 2[) . This 
increase occurs because the CO emission in NGC 2403 is 
very close to the noise limit over a large area of the galaxy. 
Finally, if these lower luminosity galaxies tend to have sys- 
tematically lower star formation rate surface densities, the 
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excitation of the CO 3-2 line relative to the ground state line 
may be lower than in the more luminous galaxies, which 
would again cause them to appear underluminous in CO 
J=3-2 (but perhaps not in CO J=l-0). 



We can use the average Lfir/L 



CO(3-2) 



ratio to make a 



crude estimate of the global molecular gas depletion time in 
the spiral galaxies in our sample. We can convert £co(3-2) 
to molecular gas mass M mo i = 1.36Mjj 2 using equation (5) 
and Lfir to sta r formation rate using equation (4) from 
iKennicuttl (|l998h . We note that this star formation rate 
equation is only strictly appropriate for starburst galaxies, 
in particular because the far-infrared luminosity may in- 
clude a co n tribu tion from dust heated by older stars; see 
IKennicuttl (|l998l ) for further details. We further convert 
from the Lir used in iKe nnicutt ( 1998 ) to Lfir by adopt- 
ing LlR = 1.3Lfir iG racia-Carpio et ail [2008) and con- 
vert to the double-p ower law initial mass function used by 
ICalzetti et all (120071 ) to give SFR = 1.3 X W~ 10 L F ir with 
SFR in Mq yr" 1 and Lfir in Lq. With these two equa- 
tions, the average Lfir/ Lco{3-2) ratio of 62±5 corresponds 
to a molecular gas depletion time of 3.0 ± 0.3 Gyr. De- 
spite the uncertainties inherent in this calculation, this mean 
value is in line with with the recent estimate of the mol ecular 
gas depletion time of 2.35 Gyr by iBigiel et~al1 l|201ll ). The 
agreement is even more striking when we consider that the 
two analyses use different CO data, different data to trace 
the star formation rate, and one is a global measurement 
while the other is based on resolved measurements. 

In cont r ast to this star formation rate analysis, 
iBendo et al.l (|2010bl . |2012| ) and others have presented ev- 
idence that much of the dust emission is at least partially 
heated by a more quiescent stellar population. Even at wave- 
lengths as short as 60 and 100 fim, there can be a sig- 
nificant contribution to dust heating from the general in- 
terstellar radiation field; one of the most spectacular ex- 
amples of this phenomenon is M31, in which 90% of the 
dust emission is not directly associate with star formation 
|Walterbos &: Schwering| |l987). Since the far-infrared lumi- 
nosity is also related to the total mass of dust, a plot of 
Lfir versus Leo may be primarily probing the molecu- 
lar gas-to-dust mass ratio in galaxies. Figure [6] shows the 
LfirI Lco('i-2) ratio plotted as a function of relative mass 
fractions of atomic and molecular gas, Mhi /Mh 2 - The two 
ratios are correlated at better than the 99% level, even when 
outliers such as M82 and NGC 4236 are included. If the far- 
infrared emission originated purely in dust associated with 
the molecular phase of the ISM traced by the CO J=3-2 line, 
we would expect no correlation in this plot. At the oppo- 
site extreme, if the far-infrared emission originated purely in 
dust associated with the atomic phase of the ISM, we would 
expect a slope of unity. The observed slope of ~ 0.2 — 0.3 
suggests that some of the far-infrared emission is originat- 
ing in dust associate with atomic gas or perhaps with more 
diffuse molecular gas that does not emit strongly in the CO 
J=3-2 line, and that this is a particularly important effect in 
galaxies where the ISM is predominantly atomic. Galaxies 
in the upper right portion of the figure tend to be either very 
early or very late-type spiral galaxies. We note that if the 
H2 mass is underestimated in some of these galaxies due to 
low metallicity, this effect would tend to steepen the slope 
and not remove the observed trend. One additional compli- 
cation in interpreting this figure is that the spatial extent 



of the atomic gas tends to be much larger than the CO or 
the far-infrared emission. The atomic gas is distributed in 
a roughly constant surface density disk, while the dust and 
molecular gas trace more of an exponential disk. 

4.3 Lco(3-2) and Lfir at low and high redshift 

llono et al l (|2009l ) presented a detailed comparison of the 
CO J=3-2 luminosity for a sample of local (D < 200 Mpc) 
luminous and ultraluminous infrared galaxies with a high- 
redshift sample of submillimeter galaxies, quasars, and Ly- 
man Break Galaxies. Including a small sample of local star 
forming galaxies, they found a strong and nearly linear cor- 
relation of Lco{3-2) an d Lfir over nearly five orders of 
magnitude, which suggested the star formation efficiency 
was constant to within a factor of two across many different 
types of galaxies and epochs. One limitation of the analy- 
sis was that the sample of local galaxies used was relatively 
small (14) and the data were limited to a single central point- 
ing (jKomu gi 2007; Mauersber gcr 1999). In addition, the far- 
infrared luminosity for the 22" diameter aperture was cal- 
culated from an H a -derived star formation rate. With the 
larger sample available here, we are in a good position to 
re-examine this relation, and to do so using global galaxy 
luminosities across all redshifts and luminosities. 



Figure [7] plots Lfir versus L 



CO(3-2) 



and also the ratio 



Lfir/L co(3 _ 2 ) for both the NGL S with the U/LIR G and 
high redshift galaxy sample from [iono et al.l |2009). This 
plot shows a significant trend of the Lfir/ Lco(3-2) ratio 
as a function of CO luminosity. The mean Lfir/ Lco(3-2) 
for galaxies with log Lfir > 11 from the local U/LIRG and 
high-redshift samples is 320 ± 40 with a standard deviation 
of 240. This is a factor of five larger than the mean ratio of 
60 measured for NGLS galaxies with 9.5 < log Lfir < 10.7. 
We can again convert the Lfir/ Lco(3-2) to a molecular 
gas depletion time, but this time using the CO-to-H2 con- 
version factor (0.8) and CO J=3-2/J=l-0 l ine ratio (0.5) 
appropriate for ULIRGs ([Wilson et al.l I2008T ) . Using equa- 
tion (3) from I Wilson et aL I J2008I) and the same star forma- 
tion rate equation from Kennicuttl (|l998l ) modified as de- 
scribed above, the Lfir/ ico(3-2) ratio of 320 corresponds 
to a molecular gas depletion time of only 50 Myr, more than 
a factor of 50 smaller than the ratio for normal disk galax- 
ies. A fit to the data for galaxies with £co(3-2) > 7.5 (or 
log Lfir > 9.5) gives a slope of ~ 0.2, which is consistent 
with a ~ 3 times higher line ratio for 10 12 L0 and a ~ 5 
times higher line ratio for 10 13 Lq than for a typical local 
disk galaxy with Lfir = 10 10 L. 



0- 



Recently, iGenzel et alj l|2010l ') have analyzed the star 
formation properties of a sample of low and high redshift 
galaxies, with one part of their analysis examining the 
Lfir/ Lco(i-o) ratio. Their study combined CO J=l-0 data 
for redshift zero galaxies with CO J=l-0 values for high red- 
shift galaxies that have been inferred from observations of 
the CO J=2-l, J=3-2, and J=4-3 transitions assuming ap- 
propriate line ratios for each type of galaxy. They also used 
a variety of star formation rate tracers, which were then 
converted to Lfir using the same equations as we have 
adopted above. A similar analysis with a some what different 
high r edshift sample has also been published by D addi et al.l 
(2010). Given the various assumptions and conversions in 



these analyses of the Lfir/L 



CO(l-O) 



ratio, it is useful to 
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Figure 6. Lfir/Lqq^_ 2 j plotted versus Mh I /Mn 2 f° r au SINGS galaxies in the NGLS with far-infrared data. Galaxies which were 
not detected in CO are omitted from this figure 



compare their results with our results which are based on 
direct measurements of both the CO and the far-infrared 
luminosities without any need for additional corrections. 

The difference in the Lfir/ 'Loo(3-2) ratio of a factor of 
5 between normal low-redshift galaxies and mergers at both 
low and high redshift is consistent with the mean difference 
of 5.4 between the average Lfir/ Lco d-o) ratio of norma l 
galaxies and merging galaxies seen bv lGenzel etafl (|20ld V 
The normal galaxy average value includes galaxies at both 
low and high redshifts. When these data are converted to a 
gas depletion time, the average for the low redshift galaxies 
is 1.5 Gyr while the average for the high-redshift galaxies 
is 0.5 Gyr. Using the average L fir/Lcq(i-q) ratio given in 
Figure 2 of lGenzel et all (|2010h . we calculate an overall aver- 
age depletion time of 1.2 Gyr for the normal galaxies across 
all redshifts, which suggests that the Lfir/ L C o(\-o) ratio of 
the low and high redshift mergers is ~ 7 times higher than 
that of normal galaxies at z=0 alone. Interesti ngly, the aver- 
age L fir/L C q(i-o) ratio given in Figure 2 of lGenzel et al.l 
(2010) for the mergers implies a gas depletion time of only 70 
MyiQ quite consistent with our estimate of 50 Myr from the 
Lfir/ Lco(3-2) ratio. Ov erall, the analysis pr esented here is 
complementary to that of lGenzel et al.l |2010f l and provides 
some evidence that the conversions between C O transitions 
and d ifferent star formation rate tracers used in lGenzel et al.l 
(|2010T ) have not biased the high-redshift analysis. 



e Note that while jGenzel et al.ll2010T) give a gas depletion time 
of 200 Myr for these data, it appears that their calculation of 
this particular number neglected to take into account the smaller 
CO-to-H2 conversion factor appropriate for merger galaxies. 



5 CONCLUSIONS 

The James Clerk Maxwell Telescope Nearby Galaxies 
Legacy Survey (NGLS) comprises an Hi-selected sample 
of 155 galaxies spanning all morphological types with dis- 
tances less than 25 Mpc. We have used new, large-area CO 
3=3-2 maps from the NGLS to examine the molecular gas 
mass fraction and the correlation between far-infrared and 
CO lum inosity using 47 galax ies drawn from the SINGS 
sample IjKennicutt et all I2003T ). We find a good correla- 
tion of th e CO J=3-2 lumi nosity with the CO J=2-l lu- 
minosity (iLerov et al l 120091 ) and the CO J=l-0 luminos- 
ity (|Kuno et all 2007 ). with average global line ratios of 
0.18±0.02 for the CO J=3-2/J=l-0 line ratio and 0.36L0.04 
for the CO J=3-2/2-l line ratio. The galaxies in our sam- 
ple span a wide range of ISM mass and molecular gas mass 
fraction, with 21% of the galaxies in the sample having more 
molecular than atomic gas and 60% having molecular to 
atomic mass ratios less than 0.5. 

We explore the correlation of the far-infrared luminos- 
ity, which traces star formation, with the CO luminosity, 
which traces the molecular gas mass. We find that the more 
luminous galaxies in our sample (with log Lfir > 9-5 or 
logL C 0(3-2) > 7.5) have a uniform Lfir/ L C o(3-2) ratio of 
62 ± 5. The lower luminosity galaxies show a wider range 
of Lfir/ Lco(3-2) ratio and many of these galaxies are not 
detected in the CO J=3-2 line. We can convert this lumi- 
nosity ratio to a molecular gas depletion time and find a 
value of 3. ± 0.3 Gyr, i n goo d agreement with recent esti- 
mates from lBigiel et all f)201 lh . We also find a correlation of 
the Lfir/ Lco(3-2) ratio with the mass ratio of atomic to 
molecular gas, which suggests that some of the far-infrared 
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Figure 7. (a) log LpjR plotted versus logL<jo(3_2) for all SINGS galaxies in the NGLS wi th far-infrared data (filled symbols and 
arrows; see Figure 0. The local U/LIRG and high redshift samples from l|lono et alJ l2009h are plotted as the open symbols, (b) 
log Lpjr/Lqq(3_2) plotted versus log £co(3— 2) f° r the same set of galaxies. 



emission originates from dust associated with the atomic 
gas. This effect is particularly important in galaxies where 
the gas is predominantly atomic. 

By comparing the NGLS data with merging gal a xies a t 
low and high redshift from the sample of llono et alJ {2009), 



which have also been observed in the CO J=3-2 line, we 
show that the Lfir/ Lco{3-2) ratio shows a significant trend 
with luminosity and thus that the correlation of Lfir with 
Lcorj -21 is not as linear a trend as found by llono et all 
(2009). Taking into account differences in the CO J=3- 
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2/J=l-0 line ratio and CO-to-Hb conversion factor between 
the mergers and the normal disk galaxies, this trend in the 
Lfir/ Lco(3-2) ratio is consistent with a molecular gas de- 
pletion time of only 50 Myr in the merger sample, roughly 
60 times shorter than in the nearby normal galaxies. 



ACKNOWLEDGMENTS 

The James Clerk Maxwell Telescope is operated by The 
Joint Astronomy Centre on behalf of the Science and 
Technology Facilities Council of the United Kingdom, the 
Netherlands Organisation for Scientific Research, and the 
National Research Council of Canada. The research of P.B., 
P. O, S. O, M.F., M.H., J.I. E.R., E.S., K.S., T.W., and 
C.D.W. is supported by grants from NSERC (Canada). 
A.U. has been supported through a Post Doctoral Re- 
search Assistantship from the UK Science & Technology 
Facilities Council. Travel support for B.E.W. and T.W. 
was supplied by the National Research Council (Canada). 
We acknowledge the usage of the HyperLeda database 
( |http: / /leda.univ- lyonl.fr} . C.V. received support from the 
ALMA-CONICYT Fund for the Development of Chilean As- 
tronomy (Project 31090013) and from the Center of Ex- 
cellence in Astrophysics and Associated Technologies (PBF 
06). This research has made use of the NASA/IPAC Ex- 
tragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration. 



REFERENCES 

Ade, P. A. R. et al., 2011, A&A, 536, 17 

Alton, P. B., Bianchi, S., Richer, J., Pierce-Price, D., & 

Combes, F., 2002, A&A, 388, 446 
Arimoto, N., Sofuie, Y., & Tsujimot, T., 1996, PASJ, 48, 

275 

Arzoumanian, D., et al., 2011, A&A, 529, L6 
Bendo G. J. et al., 2006, ApJ, 652 283 
Bendo, G. J. et al., 2007, MNRAS, 380, 1313 
Bendo, G. J. et al., 2010a, MNRAS, 402, 1409 
Bendo, G. J., et al., 2010b, A&A, 518, L65 
Bendo, G. J., et al., 2012, MNRAS, 419, 1833 
Berry, D. S., Reinhold, K., Jenness, T., & Economou, F., 
2007, in Astronomical Data Analysis Software and Sys- 
tems XVI, ASP Conference Series, Vol. 376, eds. R. A. 
Shaw, F. Hill & D. J. Bell., 425 
Bigiel, F. et al., 2008, AJ, 136, 2846 
Bigiel, F. et al. 2011, ApJ, 730, L13 

Bolatto, A., Leroy, A. K., Rosolowsky, E., Walter, F., & 

Blitz, L., 2008, ApJ, 686, 948 
Boselli, A. et al., 2010, PASP, 122, 261 
Bot C. et al., 2010, A&A, 523A, 20 

Braine, J. Combes, F., Casoli, F., Dupraz, C, Gerin, M., 
Klein, U., Wielebinski, R., & Brouillet, N., 1993, A&AS, 
97, 887 

Braine, J., et al. 1997, A&A, 326, 963 

Bregman, J. N., Snider, B. A., Grego, L., & Cox, C. V., 

1998, ApJ, 499, 670 
Buckle, J. V., et al., 2009, MNRAS, 399, 1026 



Buta, R. J., Corwin, H. G., Odewahn, S. C. 2007, The 
de Vaucouleurs Atlas of Galaxies (Cambridge: Cambridge 
University Press) 

Calzetti, D. et al. 2007, ApJ, 666, 870 

Cannon, J. M. et al., 2011, ApJ, 735, 36 

Chung, A., van Gorkom, J. H., Kenney, J. D. P., Crowl, 
H., & Vollmer, B., 2009, AJ, 138, 1741 

Combes, F. & Becquaert, J.-F., 1997, A&A, 327, 453 

Currie, M. J., Draper, P. W., Berry, D. S., Jenness, T., Ca- 
vanagh, B., & Economou, F., 2008, in Astronomical Data 
Analysis Software and Systems, ASP Conference Series, 
Vol. 394, eds. R. W. Argyle, P. S. Bunclark, & J. R. Lewis, 
650 

Dalcanton, J. J., et al., 2009, ApJS, 183, 67 

Dale, D. A., et al., 2007, ApJ, 655, 863 

Daddi, E., et al., 2010, ApJ, 714, L118 

Davies, J., et al., 2004, MNRAS, 349, 922 

Davies, J. I., et al., 2010, A&A, 518, L48 

de Vaucoulers, G., de Vaucoulers, A., Corwin, H. O, Buta, 

R. J., Paturel, G., & Fouque, P., 1991, Third Reference 

Catalogue of Bright Galaxies (New York: Springer) 
Downes, D., & Solomon, P.M. 1998, ApJ, 507, 615 
Dumke, M., Nieten, C, Thuma, G., Wielebinski, R., & 

Walsh, W., 2001, A&A, 373, 853 
Dunne, L, Eales, S., Edmunds, M., Ivison, R., Alexander, 

P., & Clements, D. L., 2000, MNRAS, 315, 115 
Eales, S. et al. 2010, A&A, 518, L62 

Elmegreen, B. G., & Elmegreen, D. M., 1986, ApJ, 311, 
554 

Foyle, K. Rix, H.-W., Walter, F., & Leroy, A. K., 2010, 

ApJ, 725, 534 
Freedman, W. L., et al., 2001, ApJ, 553, 47 
Galliano, F., Madden, S. C, Jones, A. P., Wilson, C. D., 

Bernard, J-P, & Le Peintre, F.,2003, A&A, 407, 159 
Galliano, F., Madden, S. C, Jones, A. P., Wilson, C. D., 

& Bernard, J.-R, 2005, A&A, 434, 867 
Galliano F. et al., 2011, A&A, 536, A88 
Gao, Y. & Solomon, P. M., 2004, ApJ, 606, 271 
Genzel, R. et al. 2010, MNRAS, 407, 2091 
Gil de Paz, A., et al., 2011, ApJS, 173, 185 
Giovanelli, R., et al., 2005, AJ, 130, 2613 
Gracia-Carpio, J., Garcfa-Burillo, S. Planesas, P., Fuente, 

A., & Usero, A., 2008, A&A, 479, 703 
Heifer, T. T., & Blitz, L., 1993, ApJ, 419, 86 
Heifer, T. T., Thornely, M. D., Regan, M. W., Wong, T., 

Sheth, K., Vogel, S. N., Blitz, L., & Bock, D. C.-J., 2003, 

ApJS, 145, 259 
Iono, D. et al., 2009, ApJ, 695, 1537 
Irwin, J. et al., 2011, MNRAS, 419, 1423 
Israel F.P., 1997, A&A 328, 471 
Karachentsev, I. D., et al., 2002, A&A, 383, 125 
Karachentsev, I. D., Karachentseva, V. E., Huchtmeier, W. 

K., & Makarove, D. I., 2004, AJ, 127, 2031 
Kennicutt, R. C, 1989, ApJ, 344, 685 
Kennicutt, R. C. 1998, ARA&A, 36, 189 
Kennicutt, R. C. et al. 2003, PASP, 115, 928 
Kennicutt, R. C. et al. 2007, ApJ, 671, 333 
Kennicutt, R.C., et al., 2011, PASP, 123, 1347 
Komugi, S., Kohno, K., Tosaki, T., Nakanishi, H., Onodera, 

S., Egusa, F., & Sofue, Y. 2007, PASJ, 59, 55 
Knapp, G. R., Guhathakurta, P., Kim, D.-W., & Jura, M. 

A., 1989, ApJS 70, 329 



© 0000 RAS, MNRAS 000, 000-000 



The L C o(3-2)-L FIR correlation in the SINGS sample 17 



Kuno, N. et al., 2007, PASJ, 59, 117 

Lada, E. A., Depoy, D. L., Evans, N. J., & Gatley, I., 1991, 
ApJ, 371, 171 

Lada, E. A., Bally, J., & Stark, A. A., 1991, ApJ, 368, 432 

Lee & Irwin, 1997, ApJ, 490, 247 

Lee, McCall, & Richer 2003, AJ, 125, 2975 

Leech, J., Isaak, K. G., Papadopoulos, P. P., Gao, Y., & 

Davis, G. R., 2010, MNRAS, 406, 1364 
Leonard, D. C, et al, 2002, AJ, 124, 2490 
Leroy, A. K., Walter, F., Brinks, E., Bigiel, F., de Blok, 

W. J. G., Madore, B., & Thornley, M. D., 2009, AJ, 136, 

2782 

Leroy, A. K., et al. 2009, AJ, 137, 4670 

Leroy, A. K., et al., 2011, ApJ, in press (ArXiv: 1102.4618) 

Lisenfeld U. et al., 2002, A&A, 382, 860 

Lisenfeld. U., et al., 2007, A&A, 462, 507 

Macri, L. M., Stetson, P. B., Bothun, G. D., Freedman, W. 

L., Garnavich, P. M., Jha, S., Madore, B. F. & Richmond, 

M. W., 2001, ApJ, 559, 243 
Madden, S. C, Poglitsch, A., Geis, N., Stacey, G. J., & 

Townes, C. H., 1997, ApJ, 483, 200 
Madden, S. C, Galliano, F., Jones, A. P., & Sauvage, M., 

2006, A&A, 446, 877 

Mao, R. Q., Schulz, A., Henkel, C, Mauersberger, R., Mud- 

ers, D., & Trung, D. V., 2010, ApJ, 724, 1336 
Martin, M. C, 1998, A&AS, 131, 77 

Mauersberger, R., Henkel, C, Walshi, W., & Schulz, A. 

1999, A&A, 341, 256 
Mei, S., et al. 2007, ApJ, 655, 144 

Meier, Da. S., Turner, J. L., Crosthwaite, L. P. & Beck, S. 

G, 2001, AJ, 121, 740 
Moore, B, Lake, G., & Katz, N., 1998, ApJ, 495, 139 
Mould, J. R., et al. 2000, ApJ, 529, 786 
Moustakas, J., Kennicutt, R. C, Tremonti, C. A., Dale, D. 

A., Smith, J.-D. T. & Calzetti, D., 2010, ApJS, 190, 233 
Muiioz-Mateos, J. C, et al. 2009, ApJ, 703, 1569 
Muraoka, K., et al., 2007, PASJ, 59, 43 
O'Halloran B. et al., 2010, A&A, 518, L58 
Papadopoulos, P. P., van der Werf, P., Xilouris, E. M., 

Isaak, K. G., Gao, Y. & Muehle, S., 2012, ApJ, 751, 10 
Paturel, G., Theureau, G., Bottinelli, L., Gouguenheim, L., 

Coudreau-Durand, N., Hallet, N. & Petit, C, 2003, A&A, 

412, 57 

Rosolowsky, E., Keto, E., Matsushita, S., & Willner, S. P., 

2007, ApJ, 661, 830 

Roussel, H. et al., 2010, A&A, 518, L66 
Sanders, D. B. & Mirabel, F., 1996, ARAA, 34, 749 
Sanchez-Gallego, J. R. . et al., 2011, A&A, A&A, 527, 16 
Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, J. 

A., & Soifer, B. T., 2003, AJ, 126, 1607 
Seth, A. C, Dalcanton, J. J., & de Jong, R. S., 2005, AJ, 

129, 1331 

Strong, A. W., et al. 1988, A&A, 207, 1 

Tacconi, L. J., et al. 2008, ApJ, 680, 246 

Tacconi, L. J., et al., 2010, Nature, 463, 781 

Thronson, H. A., 1988, in "Galactic and Extragalactic Star 

Formation" eds. R. E. Pudritz & M. Fich, NATO ASI 232, 

[Dordrecht: Kluwer] 621 
Tonry, J. L, Dressier, A., Blakeslee, J. P., Ajhar, E. A., 

Fletcher, A. B., Luppino, G. A., Metzger, M. R., & Moore, 

C. B., 2001, ApJ, 546, 681 
Tosaki, T., Miura, R., Sawada, T., Kuno, N., Nakanishi, 



K., Kohno, K., Okumura, S. K. & Kawabe, R., 2007, ApJ, 
664, L27 

Vlahakis, C, Dunne, L., & Eales, S. 2005, MNRAS, 364, 
1253 

Vogel, S. N., Kulkarni, S. R., & Scoville, N. Z., 1988, Na- 
ture, 334, 402 

Walter, F., Brinks, E., de Blok, W. J. G., Bigiel, F., Ken- 
nicutt, R. C, Thornley, M. D., & Leroy, A. K., 2008, AJ, 
136, 2563 

Walterbos, R. A. M., & Schwering, P. B. W., 1987, A&A, 
180, 27 

Ward, J. S., Zmuidzinas, J., Harris, A. I. & Isaak, K. G., 

2003, ApJ, 587, 171 
Warren, B. E., et al., 2010, ApJ, 714, 571 
Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 

693 

Wilson, C. D. & Walker, C. E., 1994, ApJ, 432, 148 
Wilson, C. D., 1995, ApJ, 448, L97 

Wilson, C. D., Walker, C. E., & Thornley, M. D., 1997, 

ApJ, 483, 210 1994, ApJ, 432, 148 
Wilson, C. D. et al., 2008, ApJS, 178, 189 
Wilson, C. D. et al., 2009, ApJ, 693, 1736 
Wilson, C. D. et al., 2011, MNRAS, 419, 1409 
Young, J. S. & Scoville, N. Z. 1991, ARAA, 29, 581 
Young, J. S., et al., 1995, ApJS, 98, 219 
Zavagno, A., et al., 2010, A&A, 518, L101 
Zhu M. et al., 2009, ApJ, 706, 941 



APPENDIX A: FIELD AND VIRGO GALAXY 
MEMBERSHIP IN THE NGLS 

The galaxies chosen to constitute the Field and Virgo sam- 
ples are listed in Tables IAUA3I Note that not every Virgo 
galaxy with an H I detection will be included in our sample 
because of the Hi flux limit and the constraints on the size 
of the sample. See § I2.1l for more details. 



APPENDIX B: EFFECT OF SHARED OFFS ON 
CALCULATING UNCERTAINTY FOR L co(3 _ 2) 

Both the jiggle and raster mapping modes used in this survey 
make use of a shared "off" position in order to increase the 
mapping efficiency. For the jiggle observations, the off posi- 
tion is observed after each complete cycle through the 16 dif- 
ferent pointings that produce a complete map. For a raster 
observation, the off position is observe after the completion 
of a scan of a single row. The result of using a shared off is 
that the noise is somewhat correlated between different pix- 
els in the map because the same off position is used for many 
"on" observations. While the noise is estimated properly for 
any individual pixel in the map, summing many pixels to 
obtain the total luminosity of an extended object requires 
a careful estimate of the resulting uncertainty. In effect, if 
cjj is th e measurement uncertainty in a single pixel, then 
of will underestimate the uncertainty in the summed 
luminosity because of the effect of the shared offs. 

Consider first a very simple situation where we have n 
ONs all sharing the same single OFF measurement. If we 
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Table Al. Field galaxies in the NGLS 



Name 



a(J2000.0) a 5(J2000.0) a V hel b 





(h:m:s) 




( ■ 




) 


(km s J 


ESO 538-024 


00 


:10 


:17. 


9 


-18: 


15: 


:54 


1551 


NGC 0210 


00 


:40 


:35. 


1 


-13: 


52: 


26 


1663 


NGC 0216 


00 


:41: 


:26. 


8 


-21: 


02: 


15 


1541 


IC 0051 


00 


:46: 


:24. 


3 


-13: 


26: 


32 


1753 


NGC 0274 


00 


:51 


:01. 


7 


-07: 


03: 


22 


1773 


NGC 0404 


01 


:09 


:27. 





+35: 


13: 


01 


18 


NGC 0450 


01 


:15: 


:30. 


9 


-00: 


51: 


:38 


1829 


NGC 0473 


01 


:19 


:55. 


1 


+16: 


32: 


10 


2222 


NGC 0615 


01 


:35 


:05. 


7 


-07: 


20: 


25 


1948 


ESO 477-016 


01 


:56: 


:16. 





-22: 


51: 


01 


1646 


NGC 1036 


02 


:40: 


:29. 





+ 19: 


17: 


49 


802 


NGC 1140 


02 


:54 


:33. 


6 


-10: 


01: 


:42 


1492 


NGC 1156 


02 


:59 


: 12. 


8 


+25: 


11: 


.27 


383 


ESO 481-019 


03 


:18 


: 13. 


3 


-23: 


16: 


:55 


1480 


NGC 1325 


03 


:24: 


:25. 


6 


-21: 


32: 


:38 


1646 


NGC 2146A 


06 


:23 


:54. 


3 


+78: 


31: 


50 


1386 


NGC 2742 


09 


:07 


:33. 


5 


+60: 


.28: 


15 


1285 


NGC 2787 


09 


:19: 


:19. 


1 


+69: 


12: 


12 


695 


UGC 05272 


09 


:50: 


:22. 


1 


+31: 


29: 


16 


524 


NGC 3077 


10 


:03 


:20. 


1 


+68: 


11: 


01 


-14 


NGC 3162 


10 


:13 


:31. 


6 


+22: 


11: 


15 


1384 


NGC 3227 


10 


:23: 


:30. 


7 


+19: 


51: 


51 


1131 


NGC 3254 


10 


:29: 


:20. 


1 


+29: 


29: 


32 


1262 


NGC 3353 


10 


:45 


:22. 


4 


+55: 


57: 


37 


941 


NGC 3413 


10 


:51 


:20. 


.7 


+32: 


45: 


59 


667 


NGC 3447B 


10 


:53: 


:29. 


8 


+ 16: 


17: 


02 


1023 


UGC 06029 


10 


:55: 


:02. 


3 


+49: 


13: 


33 


1403 


NGC 3507 


11 


:03 


:25. 


4 


+18: 


08: 


07 


967 


UGC 06161 


11 


:06 


: 19. 


9 


+43: 


13: 


25 


804 


ESO 570-019 


11 


:20 


:12. 


2 


-21: 


:28: 


11 


1249 


UGC 06378 


11 


:22 


:08. 


3 


+69: 


37: 


51 


1309 


UGC 06566 


11 


:35 


: 13. 


6 


+58: 


11: 


33 


1249 


NGC 3741 


11 


:36 


:05. 


8 


+45: 


17: 


03 


223 


UGC 06578 


11 


:36: 


:36. 


8 


+00: 


19: 


00 


1178 


NGC 3782 


11 


:39: 


:20. 


6 


+46: 


30: 


51 


748 


UGC 06792 


11 


:49 


:23. 


3 


+39: 


16: 


15 


832 


NGC 3931 


11 


:51 


:13. 


4 


+52: 


00: 


03 


929 


NGC 3928 


11 


:51 


: 17. 


7 


+48: 


40: 


59 


951 


NGC 3998 


11 


:57: 


:56. 


2 


+55: 


27: 


13 


1099 


NGC 4013 


11 


:58: 


:31. 


5 


+43: 


56: 


:49 


744 


IC 0750 


11 


:58: 


:52. 





+42: 


43: 


19 


694 


UGC 07009 


12 


:01 


:41. 





+62: 


19: 


40 


1097 


NGC 4041 


12 


:02 


:12. 


2 


+62: 


08: 


11 


1221 


NGC 4117 


12 


:07: 


:46. 


2 


+43: 


07: 


35 


886 


NGC 4138 


12 


:09 


:29. 


8 


+43: 


11: 


07 


908 


NGC 4190 


12 


:13 


■Al. 


2 


+36: 


:37: 


53 


192 


IC 3105 


12 


:17: 


:33. 


7 


+ 12: 


23: 


11 


-159 


NGC 4288 


12 


:20 


:38. 


3 


+46: 


17: 


31 


556 


UGC 07428 


12 


:22 


:02. 


7 


+32: 


05: 


11 


1170 


UGC 07512 


12 


:25 


Al. 


1 


+02: 


09: 


35 


1505 


NGC 4504 


12 


:32: 


:17. 


5 


-07: 


33: 


:48 


1000 


NGC 4550 


12 


:35: 


:30. 


7 


+ 12: 


13: 


15 


409 


UGC 07827 


12 


:39 


:39. 





+44: 


19: 


12 


552 


PGC 043211 


12 


:47 


:59. 


9 


+ 10: 


:58: 


32 


1140 


NGC 4772 


12 


:53: 


:29. 


1 


+02: 


10: 


05 


998 



a From RC3 (de Vaucouleurs et al. 1991) 

b Systemic velocity from the H I line (heliocentric) 
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Table A2. Field galaxies in the NGLS (continued) 



Name a(J2000.0) a 5(J2000.0) a V he i b 

(h:m:s) (°:':") (km s" 1 ) 



NGC 4941 


13 


04 


13 





-05 


33 


06 


1108 


PGC 045195 


13 


04 


31 


2 


-03 


31 


20 


1362 


UGC 08303 


13 


13 


17 


9 


+36 


12 


56 


929 


ESO 508-030 


13 


14 


55 


2 


-23 


08 


14 


1510 


NGC 5477 


14 


05 


33 


2 


+54 


27 


40 


328 


NGC 5486 


14 


07 


25 


1 


+55 


06 


10 


1368 


PGC 140287 


14 


16 


57 


3 


+03 


50 


03 


1497 


IC 1024 


14 


31 


27 


1 


+03 


00 


29 


1411 


NGC 5701 


14 


39 


11 


1 


+05 


21 


49 


1524 


IC 1066 


14 


53 


02 


9 


+03 


17 


45 


1581 


PGC 057723 


16 


17 


15 


8 


-11 


43 


54 


934 


NGC 6140 


16 


20 


58 





+65 


23 


23 


866 


NGC 6118 


16 


21 


18 


6 


-02 


17 


02 


1611 


PGC 058661 


16 


38 


08 


9 


-04 


49 


23 


1581 


IC 1254 


17 


11 


33 


5 


+72 


24 


07 


1283 


NGC 7465 


23 


02 


01 





+15 


57 


53 


1972 


NGC 7742 


23 


11 


15 


8 


+10 


46 


01 


1677 



a From RC3 (de Vaucouleurs et al. 19911 

b Systemic velocity from the H I line (heliocentric) 



average n ONs, the uncertainty in that average is given by 



0~ON,avg — 



ON,i/ r 



It is convenient to define (Nuoff) 2 ~ =1 (NjiJOFF.i) 2 
so that 



°L = ^ &ON,i + (N(TOFF) 



But since the OFF is identical for all measurements, we have 
o~oFF,avg = ooff and the uncertainty in the average differ- 
ence ON-OFF is given by 



<JON-OFF,avg 



\ i=l 



2 

OFF,j 



We can rewrite this equation as 



n(70N-OFF,avg 



' n2cr OFF,j 



(Bl) 



to arrive at an expression for ox = naoN-OFF,avg, the un- 
certainty in the sum of these n pixels. 

We now derive the appropriate formula for the uncer- 
tainty in the presence of shared offs where not every pixel 
has the same off. Suppose we are summing N measurements, 
where each measurement can be represented as the differ- 
ence between an ON measurement and an OFF measure- 
ment (modulo scaling factors such as system temperature 
and so on). So then the luminosity can be represented by 



L = 



N 

E 

i=l 



ONi-^NjOFFj 



3=1 



where there are N ONs, M OFFs, and N i = N - Let 

the uncertainty in ONi be aoN,i and the uncertainty in 
OFFj be aoFF,j- The uncertainty in L, ox, is then given 
by 



2 



N 

E 



(B2) 



3=1 



Now, we do not have access to the separate ON and 
OFF measurements, merely the difference which gives us 
the signal in each individual pixel. So we will have to be a 
little clever in getting at the appropriate uncertainty for the 
sum of many pixels. The measurement uncertainty in any 
single pixel is given by tr^ eas a = oo N ,i + °OFF,j and so 



JV JV 
' meas,z 



1 = 1 



i=l 



2 ^ ^ 2 

&ON,i + / j Nj&OFF,j 



J=l 



Now it is easy to show that 



N M 

- X " ^ Nja 2 OFF>j + (NWoff) 2 (B3) 



® L / ^meas ,i / 

i=l j = l 



Now, for a given receptor in the detector array, 
a OFF,j = afmeas = a ^ o 2 m< , as tj / Nj where a = 
l/(toff /ton + 1) and t ff and t on are the integration times 
spent on the OFF and the ON position, respectively. Simi- 
larly, aoFFj = \fa Orae-aa.j /Nj and we can write equa- 
tion (B3) as 



2 _ \ ^ 2 

®L / ^meas ,i 



i=l 



a > a 2 m eas,i + a(S a 



i=l 



me a s . i I 



and we arrive at our final equation for the uncertainty in 
the sum of N pixels in the presence of shared offs, 



2 



a ) E a ™ e 



as,i "t~ G ( 



AT 

i 



For the jiggle map mode, t a ff = 4t on and so a — 0.2. For 
the raster maps, the ratio of t ff/t on varied with the size 
and dimensions of the map; we have estimated an average 
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Table A3. Small Virgo galaxies in the NGLS 



Name 



a(J2000.0) a 
(h:m:s) 



<5(J2000.0) a V hel b 
(°:':") (kms- 1 ) 



PGC 039265 


12 


:16 


:00 


.4 


01: 


39: 


01 


2198 


NGC 4241 


12 


:18 


:00 


.1 


06: 


39: 


07 


733 


NGC 4262 


12 


:19 


:30 


.6 


1 1: 


52: 


39 


1363 


IC 3155 


12 


:19: 


: 15 


.3 


06: 


00: 


21 


2209 


NGC 4268 


12 


:19 


:47 


.2 


05: 


17: 


01 


2169 


NGC 4270 


12 


:19 


:49 


.5 


05: 


27: 


18 


2351 


NGC 4277 


12 


:20 


:03 


.8 


05: 


20: 


29 


2398 


NGC 4298 


12 


:21 


:32 


.8 


1 1: 


36: 


22 


1126 


NGC 4301 


12 


:22 


:27 


.3 


01: 


33: 


58 


1275 


NGC 4318 


12 


:22 


:43 


.4 


08: 


11: 


51 


1227 


NGC 4324 


12 


:23 


:05 


.9 


05: 


11: 


59 


1685 


NGC 4376 


12 


:25 


:18. 


.1 


05: 


: 14: 


28 


1161 


NGC 4383 


12 


:25 


:25 


.6 


16: 


28: 


12 


1641 


NGC 4390 


12 


:25 


:50 


.7 


10: 


27: 


32 


1131 


NGC 4394 


12 


:25 


:55 


.7 


18: 


12: 


50 


862 


NGC 4423 


12 


:27 


:09 


.0 


05: 


52: 


19 


1104 


IC 3365 


12 


:27: 


:11 


.6 


15: 


53: 


46 


2375 


NGC 4430 


12 


:27 


:26 


.2 


06: 


15: 


16 


1468 


UGC 07590 


12 


:28 


:18 


.8 


08: 


13: 


16 


1118 


NGC 4468 


12 


:29 


:30 


.9 


1 1: 


02: 


57 


908 


NGC 4470 


12 


:29 


:37 


.8 


07: 


■19: 


23 


2366 


NGC 4522 


12 


:33 


:39 


.7 


09: 


10: 


26 


2342 


NGC 4532 


12 


:34 


:19. 


.5 


06: 


:28: 


02 


2052 


UGC 07739 


12 


:34 


:45. 


.0 


06: 


:18: 


06 


1997 


IC 3522 


12 


:34: 


:45. 


.9 


15: 


13: 


13 


668 


NGC 4561 


12 


:36 


:08 


.2 


19: 


19: 


21 


1430 


NGC 4567 


12 


:36 


:32 


.8 


11: 


15: 


:28 


2227 


NGC 4568 


12 


:36 


:34. 


.3 


11: 


11: 


19 


2256 


IC 3583 


12 


:36: 


:43. 


,7 


13: 


15: 


32 


1039 


IC 3591 


12 


:37 


:02 


.6 


06: 


:55: 


33 


1635 


IC 3617 


12; 


:39: 


:24. 


.7 


07: 


:57: 


52 


2079 


NGC 4595 


12 


:39 


:51 


.9 


15: 


17: 


52 


604 


NGC 4639 


12 


:42 


:52 


.5 


13: 


15: 


23 


977 


NGC 4640 


12 


:42 


:57 


.8 


12: 


17: 


12 


2082 


NGC 4647 


12 


:43 


:32 


.6 


11: 


31: 


57 


1334 


NGC 4651 


12 


:43 


:42. 


.7 


16: 


23: 


36 


804 



a Prom RC3 (de Vaucouleurs et al. 1991) 

6 Systemic velocity from the H I line (heliocentric) 



value a = 1.07 for the rasters and have used this in the noise 
calculations for all the galaxies. 



Hon]) using the same processing techniques used in this 
paper and so are not reproduced here. 



APPENDIX C: IMAGES OF THE SINGS 
GALAXIES FROM THE NGLS 

In this appendix, we present the images derived from the 
CO J=3-2 observations of the NGLS galaxies that a re also 
members of the SINGS sample (|Kennicutt et al.ll2003l ). Only 
galaxies with significant detections are shown here; upper 
limits to the CO luminosity for the remaining galaxies are 
given in Table [3] 

Two images are shown for each galaxy: the CO J=3-2 
integrated intensity image, and the CO contours overlaid 
on an optical image from the Digitized Sky Survey. For 
those galaxies with sufficient signal-to-noise, we also show 
the first moment map derived from the CO data cube, which 
traces the velocity field and the second moment map, which 
traces the velocity dispersion. The second moment maps for 
nine galaxies with good signal to noise and in clinations less 
than 60° have been presented and discussed in I Wilson et al.l 
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Figure CI. CO J=3-2 images for NGC 0628. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2) K km s _1 (T MS ). 
(b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity field as traced by 
the CO J=3-2 first moment map. Contour levels are (624, 632, 640, 648, 656, 664, 672, 68 0) km s" 1 . The thick line shows the systemic 
velocity given in Tabled] Si milar images derived from the same data have been published in I Warren et al The velocity dispersion 

map has been published in lWilson et al. 




30 25 2:27:20 15 10 05 30 25 2:27:20 15 10 05 

Right ascension (hh:mm:ss.s) Right ascension (hh:mm:ss.s) 

Figure C2. CO J=3-2 images for NGC 0925. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, f, 2, 4) K km s _1 
(Tms). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (525, 540, 555, 570, 585) km s _1 . 
(d) The velocity dispersion a v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16) km s . 
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Figure C3. CO J=3-2 images for NGC 2403. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4) K km s _1 
(Tjib). (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity field as 
traced by the CO J=3-2 fi rst moment map. Co ntour levels are (21, 46, 71, 96, 121, 146, 171, 196, 221) km s -1 . The velocity dispersion 
map has been published in lWilson et al.l j201ll ). 
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Figure C4. CO J=3-2 images for NGC 2976. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4) K km s _1 
(T MB ). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (-50, -32, -14, 4, 22, 40, 58, 76, 
94) km s — 1 . (d) The velocity dispersion a v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8) km s — 1 . 
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Figure C5. CO J=3-2 images for NGC 3031. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4) K km s _1 
(T M b). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (-328, -258, -188, -118, -48, 22, 92, 
162) km s — 1 . (d) The velocity dispersion cr v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16, 32, 64, 128) 
km s ■ Images derived from the same data using careful flagging and analysis to detect weak but real features have been published in 
ISanchez-Gallego et all i201ll '). 
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Figure C6. CO J=3-2 images for NGC 3034. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32, 64, 
128, 256, 512) K km s" 1 (T MS ). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (129, 
154, 179, 204, 229, 254, 279) km s" 1 . (d) The velocity dispersion cr v as traced by the CO J=3-2 second moment map. Contour levels 
are (4, 8, 16, 32, 64) km s _1 . 
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Figure C7. CO J=3-2 images for NGC 3184. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4) K km s _1 
(Tjvfs). (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity field as 
traced by the CO J=3-2 fi rst moment map. Co ntour levels are (532, 546, 560, 574, 588, 602, 616, 630) km s — 1 . The velocity dispersion 
map has been published in lWilson et al.l j201ll ). 
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Figure C8. CO J=3-2 images for NGC 3198. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8) K km s" 1 
(Ta^s). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (606, 636, 666, 696, 726, 756) km 
s — . (d) The velocity dispersion cr v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16, 32) km s — 1 . 
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Figure C9. CO J=3-2 images for NGC 3351. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32) K 
km s" 1 (Tjffl). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (702, 735, 768, 801, 834, 
867, 900) km s -1 . (d) The velocity dispersion a v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16, 32, 64) 
km s _1 . 
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Figure C10. CO J=3-2 images for NGC 3521. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8) K km s _1 
(T MB ). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (558, 613, 668, 723, 778, 833, 888, 
943, 998) km s _1 . (d) The velocity dispersion u v as traced by the CO J=3-2 second momen t map. Contour levels are (4, 8, 16, 32) km 
s — 1 . Similar images derived from the same data have been published in lWarren et al 
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Figure Cll. CO J=3-2 images for NGC 3627. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32) K 
km s — 1 (Tmb)- (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (559, 605, 651, 697, 743, 
789, 835, 881) km s — . (d) The velocity dispersion a v as traced by the CO J=3-2 se cond moment map. C ontour levels are (4, 8, 16, 32, 
64, 128) km s _1 . Similar images derived from the same data have been published in lWarren et al 
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Figure C12. CO J=3-2 images for NGC 3938. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2) K km s _1 
(Tjks). (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity field as 
traced by the CO J= 3-2 first moment ma p. Contour levels are (769, 781, 793, 805, 817, 829, 841) km s — 1 . The velocity dispersion map 
has been published in I Wilson et al. | j201ll) . 
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Figure C13. CO J=3-2 images for NGC 4236. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2) K km s 
(Tmb)- (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (-190, -150, -110, -70, -30, 10, 50, 
90, 130, 170) km s _1 . (d) The velocity dispersion <r„ as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16, 32, 
64) km s _1 . 
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Figure C14. CO J=3-2 images for NGC 4254. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16) K km 
s _1 (Tjvfs). (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity field as 
traced by the CO J=3-2 first moment map. Contour levels are (2304, 2340, 2376, 2412, 2448, 2484) km s _1 . Si milar images d e rived from 
the same data have been published in lWilson et all J200ST) ■ The velocity dispersion map has been published in I Wilson et al" 




Figure C15. CO J=3-2 images for NGC 4321. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32) K 
km s (Tmb)- (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity 
field as traced by the CO J=3-2 first moment map. Contour levels are (1 424, 1459, 1494, 152 9, 1564, 1599, 1634, 1669, 1704, 1739) 
km s -1 . Si milar images derived from the same data have been published in Wils on et al. The velocity dispersion map has been 

published in I Wilson et al.l ll201ll 1. 




Figure C16. CO J=3-2 images for NGC 4450. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1) K km s _1 (Tmb). 
(b) CO J=3-2 overlaid on a Digitized Sky Survey image. 
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Figure C17. CO J=3-2 images for NGC 4559. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2) K km s _1 
(Tats), (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (735, 758, 781, 804, 827, 850, 873, 
896) km s — 1 . (d) The velocity dispersion a v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16) km s — . 
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Figure C18. CO J=3-2 images for NGC 4569. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16) K km 
s -1 (Tms). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (-379, -339, -299, -259, -219, 
-179, -139, -99, -59, -19) km s —1 . (d) The velocity dispersion a v as traced by the CO J=3-2 second mome nt map. Contour levels are (4, 
8, 16, 32) km s _1 . Similar images derived from the same data have been published in lWilson et al.l l|2009h . 
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Figure C19. CO J=3-2 images for NGC 4579. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4) K km s~ x 
(T M b). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (1340, 1390, 1440, 1490, 1540, 
1590, 1640, 1690) km s" 1 . (d) The velocity dispersion <r„ as traced by the CO J=3-2 second moment m ap. Contour levels ar e (4, 8, 
16, 32) km s" 1 . Images of this galaxy derived from early observations in the survey have been published in lWilson et al ] (2009). These 
images are made using later data which had a higher quality. 
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Figure C20. CO J=3-2 images for NGC 4631. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32) K 
km s — 1 (Tmb)- (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (520, 550, 580, 610, 640, 
670, 700, 730, 760) km s — ■ (d) The velocity dispersion a v as traced by the CO J=3-2 sec ond momen t map . Contour levels are (4, 8, 
16, 32, 64, 128) km s _1 . Similar images derived from the same data have been published in llrwin et al ] j201l|y 




Figure C21. CO J=3-2 images for NGC 4736. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32) K 
km s — 1 (Tmb)- (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity 
field as traced by the CO J=3-2 first moment map. Con tour levels are (214, 241, 268, 295, 322, 349, 376, 403) km s" 1 . The velocity 
dispersion map has been published in Wilson ct al. ( 20111 ), 
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Figure C22. CO J=3-2 images for NGC 4826. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32, 64) 
K km s — 1 (Tmb). (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity 
field as traced by the CO J=3-2 first moment map Conto ur levels are (230, 270, 310, 350, 390, 430, 470, 510, 550) km s _1 . The velocity 
dispersion map has been published in I Wilson et al. | feOllI ). 
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Figure C23. CO J=3-2 images for NGC 5033. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16) K km 
s_1 (Tmb)- (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (593, 663, 733, 803, 873, 943, 
1013, 1083) km s~ 1 . (d) The velocity dispersion a v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16, 32, 
64, 128) km s _1 . 
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Figure C24. CO J=3-2 images for NGC 5055. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16) K km 
s _1 (Tms)- (b) CO J=3-2 integrated intensity contours overlaid on an optical image from the Digitized Sky Survey, (c) Velocity field 
as traced by the CO J=3-2 first mo ment map. Contour levels are (315, 360, 405, 450, 495, 540, 585, 630, 675) km s _1 . The velocity 
dispersion map has been published in I Wilson et alj ll201l[ ). 
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Figure C25. CO J=3-2 images for NGC 5194. (a) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4, 8, 16, 32) K 
km s~ J (T MS ). (b) CO J=3-2 overlaid on a Digitized Sky Survey image, (c) Velocity field. Contour levels are (392, 412, 432, 452, 472, 
492, 512, 532) km s -1 . (d) The velocity dispersion a v as traced by the CO J=3-2 second moment map. Contour levels are (4, 8, 16, 32, 
64) km s" 1 . A more complete analysis fo the data for NGC 5194 will be published in Vlahakis et al., in preparation. 
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Figure C26. CO J=3-2 images for the small SINGS galaxies detected in jiggle map mode. See Figure ICll for more details, (left of 
pair) CO J=3-2 integrated intensity image. Contours levels are (0.5, 1, 2, 4) K km s _1 (Tmb). (right of pair) CO J=3-2 overlaid on a 
Digitized Sky Survey image. 
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